PLASMA  ACCELERATION  IN  PULSED  ABLATIVE  ARC 

DISCHARGES 


Final  Report  For  Period  15  March  1974  - 31  March  1977 


AUTHORS:  DOMINIC  J.  PALUMBO 

MARTIN  BEGUN 


10ER  CONTRACT  NO.  Fl»  1*620- 7^«-C -0055 


APPROVED  FOR  PUBLIC  RELEASE 
DISTRIBUTION  UNLIMITED 


PREPARED  FOR 


OF  SCIENTIFIC  RESEARCH  (AFSC) 


wOSR  - TRm  77  r 0628 


PREPARED  BY: 

FAIRCHILD  INDUSTRIES,  INC. 

FAIRCHILD  REPUBLIC  COMPANY 
FARMINGDALE.N.Y.  N.Y.  11735 


NOTICE 


Reproduction,  translating  publication,  use  and  disposal  tn  whole  or 


FOREWORD 


c.,  Fairchild  Republic 
Plasma  Acceleration  tn 


Compaq?,  under  Air  Force  Contract  F44620-74  -C-0055 


MS147V003 


Research,  This  program  was  monitored  by  Dr,  Bernard  T.  Wolfson 


William  J, 


F=X\URGMILD 


CONTENTS 


Section  Page 

1.0  INTRODUCTION  1 

2.  0 PHYSICAL  BACKGROUND  AND  EXPERIMENTAL 

APPARATUS  2 

2.1  General  Description  of  Pulsed  Ablative  Plasma  Accelerators  2 

2. 2 Description  of  Laboratory  Test  Facility  3 

3. 0 DIAGNOSTICS  6 

3. 1 Total  Discharge  Current  and  Voltage  6 

3. 2 Results  of  High  Speed  Photography  7 

3. 3 Arc  Velocity  and  Shape  Using  Time-of-Arrival  Probe  13 

3.4  Distribution  of  Magnetic  Flux  Density  15 

3. 5 Current  Distribution  23 

3.6  Electron  Number  Density  and  Temperature  29 

3.6.1  Description  of  Langmuir  Probe  29 

3.6.2  Data  Acquisition  Procedure  31 

3.6.3  Data  Reduction  32 

4 3.6.4  Number  Density  and  Temperature  Distribution  34 


f 


« 


4.0  THEORETICAL  APPROACH  38 

4.1  Physical  Model  38 

4. 1. 1 Initial  Conditions  38 

4.1.2  Boundary  Conditions  40 

4.1.3  Plasma  Coupling  to  Discharge  Circuitry  46 

4.2  Governing  Equations  48 

4.3  Chemical  Model  49 

4.4  Numerical  Integration  Scheme  53 

5.0  RESULTS  OF  CALCULATIONS  55 

6. 0 CONCLUSIONS  63 

7.0  REFERENCES  65 

APPENDIX  A 67 

APPENDIX  B 70 


1 iccmm  to 

lilt 

am*  mm 

tf 

ice 

Ml  MM 

o 

iMrxo’jccn 

□ 

‘ ;VS!'HGATIM. 


litx/AHiuimn  e##o 


/ 


!L  md/Bf  SPECIAL 


! 


D D C 

ir?ror?nn  nr?f 

nay  io  ran 

tEMsnnsI 

D 


ii 


I 

II 

i 

j 

j 


OFFICE  OF  SCIENTIFIC  BESEECH  <»■» 

transmittal  to  ddc  revie,ed  and  13 

rx-*1  ia*  418 190-12 

ioa  is  unli“lted* 


Information  QfficeT 


FXXURCHILD 

' >(  **LJHL. >C  C(JM»4/VV 


ILLUSTRATIONS 


Schematic  Cross-Section  of  Pulsed  Ablative  Accelerator 
Test  Facility 

Discharge  Current,  Voltage,  and  Energy  Input  as  a 
Function  of  Time 

Mirror  Images  of  Front  and  Side  Views  of  Accelerator 

High  Speed  Photographs  of  Accelerating  Plasma 

Arc  Front  Location  and  Shape  Obtained  Using  Time-of- 
Arrival  Probe 

Small  Coil  Magnetic  Probe  Magnified  20X 

Magnetic  Probe  Response  and  Discharge  Current  Versus 
Time  at  Various  Positions  Along  the  Ablator  Surface 

Magnetic  Flux  Density  Distribution  at  Various  Times 
During  Discharge 

Current  Distribution  at  Various  Times  During  Discharge 

Langmuir  Probe  Cbnfiguration 

Schematic  of  Probe  Cleaning  Circuitry 

Langmuir  Probe  Equipped  for  Remote  Cleaning  by 
Electron  Bombardment 

Typical  Volt-Ampere  Langmuir  Probe  Characteristic 

Electron  Number  Density  as  a Function  of  Distance  Along 
the  Accelerator  Centerline  at  Three  Different  Times 

Typical  Langmuir  Probe  Response  and  Discharge  Current 
2cm  Past  the  Electrode  Exit  Plane 

Electron  Number  Density  Across  the  Interelectrode  Gap 
a Distance  2cm  Past  the  Electrode  Exit  Plane 

Equivalent  Circuit  Schematic 

Typical  Wave  Propagation  Diagram 

Computed  Magnetic  Flux  and  Current  Density  During  the 
First  Quarter  Cycle 


OOMMNV 


1. 0 INTRODUCTION 

Spacecraft  and  associated  communications,  photographic,  and  computer  technology 
4 have  reached  a level  of  sophistication  requiring  attitude  control  systems  capable  of 

maintaining  fractions  of  an  arc  second  pointing  accuracy.  Among  the  candidates 
♦ for  the  auxiliary  propulsion  subsystem  are  hydrazine,  cold  gas,  electrostatic  and 

plasma  thrusters.  Of  these,  only  the  pulsed  ablative  plasma  accelerator  offers 
highly  repeatable  discrete  impulse  bits  of  microsecond  duration  as  well  as  equivalent 
steady  state  thrust  capability.  In  addition  to  this  feature,  the  intrinsic  simplicity  and 
corresponding  reliability  of  this  concept  have  inspired  satellite  manufacturers  to 
develop  flight  qualified  microthruster  hardware,  1,2‘3’4  and  government  agencies 
to  continue  development  of  a millipound  thruster  system.  5* 6*  7 Although  this 
current  area  of  development  is  the  primary  one,  other  applications  for  plasma 
accelerators,  i.  e. ; isotope  separators,  simulators  and  injectors,  have  all  prompted 
this  basic  study  into  the  dynamics  of  the  plasma  acceleration  process  in  pulsed 
ablative  discharges. 

• To  the  present  day,  no  truly  basic  understanding  of  the  processes  of  propellant 

ablation,  plasma  production,  energy  transfer  to,  and  acceleration  of  plasma  exists  for 
this  device.  Sufficient  empirical  data  correlations  which  have  been  generated 
in  past  years  have  enabled  design  of  thrusters  to  meet  specific  requirements,  8 but 
none  of  these  include  optimization  of  thruster  efficiency.  It  is  believed  that  if  a good 
understanding  of  the  basic  processes  were  available,  such  optimization  could  be 
accomplished.  Furthermore,  other  associated  application  oriented  details,  such 
as  exhaust  plume  characteristics,  could  be  extrapolated  and,  perhaps,  predicted. 

Previous  investigations  into  pulsed  acceleration  of  plasma  have  all  been  concerned 
with  accelerating  a gaseous  working  fluid.  At  the  onset  of  this  program  it  was 
presumed  that  the  basic  processes  associated  with  accelerating  a plasma  produced 
from  a subliming  ablator  were  considerably  different  than  those  in  the  case  of  gaseous 
injection  or  static  breakdown  induced  discharges.  This  supposition  was  justified 
during  the  first  phase  of  this  effort, 9 wherein  a one-dimensional  continuum  model 
was  formulated  using  the  unsteady  equations  of  magnetohydrodynamics  as  a basis. 
Obtaining  a numerical  solution  to  those  equations  was  plagued  by  stability  problems 
and  a peculiar  non  uniqueness  associated  with  the  initial  time  step.  Nevertheless, 
the  experimental  data  and  some  analytic  work  with  MHD  equations  had  led  to  a 
clearer  picture  of  the  acceleration  mechanism. 
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During  the  second  phase  of  this  effort, 10  measurements  made  in  the  first  phase 
were  extended  to  give  a two-dimensional  picture  of  the  plasma  as  a function  of  time. 
Magnetic  flux  density,  local  average  current  density,  and  time  of  arrival  probes 
were  utilized  in  order  to  measure  some  of  the  more  salient  plasma  parameters. 

High  speed  photography  (i.e. , 2M  fps)  was  also  used  to  determine  the  motion  of  the 
luminous  plasma  front.  The  diagnostic  study  was  augmented  by  a two-species  con- 
tinuum MHD  model  in  which  the  one-dimensional  model  of  the  first  phase  was 
extended  to  two— dimensions  and  the  effects  of  finite  viscosity  and  heat  conductivity 
were  included.  At  the  end  of  that  phase  computational  difficulties  were  still  being 
encountered  but  the  reasons  for  occurrence  of  instabilities  in  the  numerics  were 
becoming  apparent  and  steps  were  being  taken  to  circumvent  these  instabilities. 

The  two-dimensional  diagnostics  were  very  successful  and  yielded  further  qualitative 
understanding  of  the  ablative  plasma  accelerator. 


The  third  and  final  phase  of  this  program  was  geared  toward  obtaining  a better  under- 
standing of  the  chemistry  involved  in  producing  plasma  by  depoiymerization  of  solid 
tetrafluoroethylene  (Teflon,  as  produced  by  DuPont).  This  study  included  Langmuir 
probe  measurements  of  the  two-dimensional  distribution  of  electron  number  density 
and  temperature,  a detailed  evaluation  of  chemical  and  ionization  reactions  possible 
among  the  products  of  depolymerized  Teflon,  and  incorporation  of  a finate  rate 
ionization  reaction  and  neutral  species  in  the  analytic  model.  The  results  of  this 
final  phase  and  of  the  previous  two  phases  are  included  in  this  final  report. 


2.0  PHYSICAL  BACKGROUND  AND  EXPERIMENTAL  APPARATUS 

2. 1 General  Description  of  Pulsed  Ablative  Plasma  Accelerators 

Energy  storage  capacitors  are  connected  in  parallel  through  a low  inductance  strip 
line  assembly  to  two  parallel  planar  electrodes  (see  Figure  1). 
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Figure  1.  Schematic  of  Pulsed  Plasma  Accelerator 
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A solid  propellant  (Teflon)  is  positioned  between  the  electrodes  at  the  same  end  at 
which  the  strip  lines  are  connected.  The  cathode  is  fitted  with  an  "igniter  plug" 
which  consists  of  a cylindrical  anode  and  cathode  separated  by  a semiconducting 
material.  The  dimensions  of  this  plug  are  much  smaller  than  those  of  the 
accelerator  electrodes  and  the  cathode  of  the  plug  is  connected  to  the  cathode  of 
the  accelerator  through  one  ohm  of  resistance.  The  accelerator  is  fired  by  first 
charging  the  capacitor  bank  to  some  voltage  VQ  which  will  appear  simultaneously 
across  the  electrodes.  The  presence  of  vacuum  between  the  electrodes,  however, 
allows  retention  of  this  charge  without  breakdown  with  the  usual  voltages  and 
interelectrode  spacings  used.  The  capacitors  are  then  discharged  through  the 
electrodes  by  firing  the  ignitor  plug.  This  is  accomplished  by  discharging  a separ- 
ate 2 microfarad  capacitor  from  600  VDC  through  a 1:3  step-up  transformer  by 
triggering  an  SCR  in  series  with  this  capacitor  on  the  primary  side.  The  ignitor 
plug,  when  fired,  produces  a miniscule  amount  of  ionized  matter  having  sufficient 
conductivity  to  break  down  the  accelerator  electrode  gap.  Since  the  one  ohm 
resistance  represents  a high  impedance  path,  the  main  capacitor  discharge  initiates 
between  the  accelerator  anode  and  cathode,  rather  than  returning  to  ground  through  the 
ignitor  plug  cathode.  Some  time  later  (i.  e. , appromimately  50-200  ns)  the  Teflon  begins 
to  depolymerize,  a process  which  transforms  the  solid  directly  into  a gas  composed 
of  heavy  molecules.  Subsequent  dissociation  and  ionization  of  this  gas  produces  a 
plasma,  into  which  the  energy  stored  in  the  capacitor  bank  is  deposited  through 
Joule  heating  and  the  creation  of  magnetic  field. 

2.2  Description  of  Laboratory  Test  Facility 

The  accelerator  and  vacuum  pumping  station  were  designed  to  allow  maximum 
flexibility  and  instrumentation  ease.  The  vacuum  station  consists  of  a glass  T-jar 
placed  on  top  of  a 4-inch  (10. 16  cm)  diameter  pump  system.  Ultimate  pressure  in 
this  facility  is  roughly  2 x 10  ^ mm  Hg.  The  main  discharge  capacitor  bank  and 
discharge  initiating  circuitry  were  located  outside  of  the  jar  to  yield  a maximum 
amount  of  space  inside.  The  accelerator  electrode  assembly  was  connected  to  the 
capacitor  bank  through  two  coaxial  feed-throughs  fabricated  from  oxygen-free 
copper.  Electrodes  measure  4. 0 inches  (10. 16  cm)  in  width  and  are  4. 0 inches 
long.  The  interelectrode  gap  is  3.0  inches  (7.62  cm)  and  the  electrodes  are  parallel 
to  each  other  over  their  length.  The  Teflon  propellant  consists  of  a 0. 313  inch 
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(7.  95  mm)  thick  slab  having  a height  of  2.  9 inches  (7.  37  cm)  anti  width  or  1. 0 inch 
(2.  54  cm).  This  slab  is  positioned  at  the  center  of  the  electrode  width  at  the  back 
end.  Two  pieces  of  0.25  inch  (6.35  mm)  Mykroy  1100  ceramic  are  placed  on  either 
side  of  the  propellant  slab  to  shield  the  remaining  electrode  area.  The  excessive 
electrode  width  (as  compared  to  propellant  width)  was  incorporated  to  eliminate 
fringe  effects  of  the  electromagnetic  fields  between  the  electrodes  in  the  vicinity 
of  the  plasma,  thereby  making  the  assumption  of  two  dimensionality  implied 
throughout  this  study  more  consistent. 

The  ignitor  plug  is  located  in  an  insulated  housing  immediately  above  the  cathode. 

The  plug  fires  into  the  interelectrode  gap  through  a 0. 25  inch  (6. 35  mm)  diameter 
chamfered  hole  in  the  cathode.  A photograph  of  the  test  facility  is  included  in 
Figure  2. 

In  order  to  facilitate  two  dimensional  scanning  of  the  plasma,  all  probes  which 
were  used  in  this  study  were  fastened  to  an  assembly  which  could  be  moved  from 
outside  the  vacuum  station  in  both  the  axial  and  latteral  directions.  This  was 
accomplished  using  rotary  motion  vacuum  feed  throughs,  flexible  shafting  and  two 
Unislide  assemblies.  This  arrangement  allowed  a complete  scan  of  the  interelectrode 
gap  to  be  made  with  each  probe  without  ever  having  to  break  the  vacuum  for 
repositioning.  It  also  allowed  extremely  accurate  positioning  of  the  probes  relative 
to  a fixed  known  origin. 

A high  degree  of  resolution  was  desired  so  as  to  determine  the  space/time  evolution 
of  the  measured  plasma  parameters  as  accurately  as  possible.  This  requirement 
implies  that  the  interelectrode  gap  be  scanned  over  a relatively  fine  grid  which 
means  that  many  individual  pieces  of  data  would  be  obtained.  Use  of  a normal 
analog  oscilloscope  would  require  taking  photographs  of  the  events  at  each  probe 
location  and  later  reducing  this  data  by  measuring,  as  accurately  as  possible,  the 
amplitude  of  recorded  signal  and  its  time  relative  to  some  consistent  reference  time 
during  the  discharge.  This  procedure  would  have  been  expensive  considering  the 
cost  of  film  and  would  also  have  required  an  excessive  man  hour  effort  for  the 
reduction  of  data  because  of  the  large  amount  of  data  to  be  accumulated.  For  this 
reason  the  decision  was  made  at  a point  during  the  program  to  rent  a digital  memory 
oscilloscope.  The  advantages  of  this  kind  of  oscilloscope  over  the  analog  type  are 
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two  fold.  Data  is  received  and  automatically  converted  to  numerical  values  in 
volts  via  an  A to  D converter.  The  analog  signal  is  sampled  at  a preset  interval  of 
time  between  samplings  and  the  data  is  stored  in  numerical  form  within  the 
built-in  memory.  The  data  is  also  displayed  in  the  form  of  dots  on  a screen  which 
are  separated  by  the  preset  time  interval,  and  this  display  remains  on  the  screen  until 
it  is  deliberately  erased  by  a new  signal  or  by  the  user.  To  retrieve  the  digital 
value  of  the  signal  at  any  of  these  points  it  is  merely  necessary  to  call  out  the  point 
in  question  using  a cursor  which  moves  from  one  point  to  the  next  and  read  off  its 
value  in  volts  (accurate  to  within  ± 1%)  directly  from  the  digital  display  which  is 
also  on  the  screen.  The  two  advantages  are,  therefore,  that  no  photographic 
record  of  the  event  is  required  since  the  data  is  obtained  directly  from  the  scope 
and  no  lengthy  expendure  of  time  is  necessary  to  generate  the  numerical  data  from  an 
analog  record.  The  savings  in  cost  of  film  alone  defrayed  approximately  half  the 
cost  of  renting  the  scope  which  is  manufactured  by  NICOLET  Instrument  Corp. 

(Model  1090A  with  93A  plug-in).  The  digital  memory  oscilloscope  has  a frequency 
response  of  1 MHz  which  is  sufficiently  high  to  record  the  events  measured  during 
this  program  since  the  discharge  cycle  has  a ringing  frequency  of  roughly  16  kHz. 
Other  plug-in  units  are  available  which  have  10  MHz  frequency  response. 

Use  of  the  digital  scope  allowed  scanning  over  a much  finer  grid  than  was  originally 
anticipated  because  of  the  large  savings  in  data  reduction  time,  and  has  contributed 
to  the  overall  success  of  the  experimental  effort  and  accuracy  of  the  data  presented 
herein. 

3. 0 DIAGNOSTICS 

3.1  Total  Discharge  Current  and  Voltage 

The  discharge  current  and  voltage  across  the  electrode  gap  were  measured  as  a 
function  of  time.  Current  was  measured  using  a Rogowski  Coil  wound  so  as  to  fit 
around  the  cathode  behind  the  propellant.  The  coil  was  calibrated  using  a precali- 
brated coil  of  known  accuracy.  Discharge  voltage  was  measured  using  a Tektronix 
1000:1  voltage  probe.  Initial  voltage  on  the  capacitor  bank  was  set  precisely  using  a 
digital  DC  voltmeter  having  an  accuracy  of  ± .1%.  The  current  and  voltage  wave- 
forms were  displayed  simultaneously  on  a Tektronix  Model  555  dual  beam  oscillo- 


scope. 


F^AlFtCHIL-D 

• 4!  t;fULK  « <>K*>-A*4Y 


The  oscillograms  are  reproduced  in  Figure  3 where  discharge  current  and  voltage 
throughout  the  first  13  microseconds  of  the  discharge  are  plotted  as  a function  of 
time.  The  energy  delivered  by  the  capacitor  bank,  given  as  the  integral  of  the 
instantaneous  power  (current  times  voltage),  is  also  plotted. 

The  current  and  voltage  curves  deviate  considerably  from  what  one  would  expect  to 
see  in  an  LRC  circuit  with  fixed  discrete  elements.  The  deviations  are  due  to  coupling 
of  the  plasma  motion  via  the  induced  magnetic  fields  and  the  time  variant  plasma 
resistance  and  inductance.  The  degree  of  deviation  suggests  strong  coupling;  i.  e. , 
the  plasma  resistance  and  inductance  are  much  larger  than  those  of  the  fixed  portion 
of  the  discharge  path.  The  fact  that  virtually  aU  of  the  energy  is  delivered  by  the 
capacitor  during  roughly  the  first  quarter  cycle  is  of  particular  interest.  Also  note 
that  subsequent  cycles  only  transfer  from  2 to  2. 5J  back  and  forth  to  the  capacitors. 

The  significance  of  this  observation  becomes  clearer  as  more  information  appears 
later  in  this  report. 

Fairly  good  energy  delivery  efficiency  (93. 1%)  from  the  capacitor  bank  indicates  a 
time-averaged  ESR  of  approximately  2.9  milliohms.  The  initial  inductance  at  the 
instant  breakdown  occurs  can  be  estimated  from  the  relationship 

-cHA  _ Lo 
' dt  t = O Lo 

Where  I is  the  total  current,  t is  time,  VQ  is  the  initial  voltage,  and  Lq  is  the  initial 
inductance  (i.  e. , of  the  circuit  plus  arc).  The  oscillograms  indicate  a dl/dt  of 
2. 6 x 1010  A/s,  yielding  a value  of  115. 6 nH  for  LQ. 

3.2  Results  of  High  Speed  Photography 

A Beckman  & Whitley  Model  189  High  Speed  Framing  Camera  was  used  to  photograph 
the  discharge  from  its  initiation  to  near  completion.  This  camera  will  take  a maximum 
of  25  exposures  at  a maximum  rate  of  one  exposure  every  0.4  microseconds.  Since  it 
was  desired  to  photograph  at  least  the  first  cycle  of  the  discharge  (i.e. , roughly  7 
microseconds)  with  a meaningful  degree  of  resolution,  a framing  rate  of  one  exposure 
per  0. 5 microseconds  was  chosen. 
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The  arc  was  simultaneously  photographed  from  the  side  and  front  using  a system  of 
four  front  surface  mirrors  to  project  both  images  side  by  side  at  the  same  level. 

A photograph  of  this  arrangement  is  presented  in  Figure  4.  The  camera  lens  was 
located  some  25  feet  (7.62m)  from  the  mirrors.  A front  view  of  the  arc  evolution  was 
desired  so  that  the  relative  lateral  velocity  of  the  plasma  could  be  determined  for  com- 
parison to  its  forward  progression.  With  this  information  the  assumption  of  two  dimen- 
sional motion  could  be  cross-checked. 

Considerable  effort  was  required  in  order  to  initiate  framing  at  the  instant  of  break- 
down. To  do  this  it  was  necessary  to  fire  the  thruster  using  the  output  signal  generated 
by  the  camera.  In  addition,  the  thruster  had  to  be  completely  isolated  from  ground 
subsequent  to  charging  the  capacitor  bank.  Several  runs  were  made  in  which  the 
angular  location  of  top-dead-center  in  the  camera  was  adjusted  until  only  a single  frame 
appeared  prior  to  initiation. 

Photography  was  done  without  the  use  of  filters  ahead  of  the  lens.  A typical  sequence 
is  presented  in  Figure  5.  Photographs  of  the  front  view  indicate  that  the  arc  begins 
as  a fairly  straight  line  centered  on  the  electrodes  and  against  the  face  of  the  Tcf*  n 
propellant.  At  0.5  us  the  line  appears  to  be  approximately  1.7  cm  in  width  and  iform. 
A half  microsecond  later  the  line  is  roughly  3. 2 cm  wide  and  at  1.5  us  it  is  3.4  cm 
wide.  The  width  appears  to  remain  constant  past  that  point  in  time. 

The  view  from  the  side  shows  that  the  arc  remains  against  the  propellant  surface 
as  the  luminous  front  moves  forward  up  until  about  2.5  us  where  a thin  dark  region  is 
visible  between  the  initial  luminous  plasma  and  a second  luminous  plasma  present  at 
the  propellant  surface.  Note  that  at  this  time  the  first  quarter  cycle  of  the  discharge 
current  has  not  yet  ended,  but  peak  current  has  been  reached  and  the  current  is  now 
decreasing.  As  time  progresses,  the  dark  region  increases  in  size.  The  initial  plasma 
continues  its  forward  motion  and  the  second  luminous  region  continues  to  remain 
against  the  propellant  surface  and  becomes  slightly  thicker.  At  zero  current  (roughly 
3.9  us)  there  is  no  discernible  irregularity  introduced  and  the  above  events  proceed 
continuously  up  until  the  exposure  at  5 us,  where  the  second  luminous  region  appears 
to  have  moved  off  the  ablator  surface  and  move  forward  behind  the  first.  A dark 
region  is  again  evident  between  this  second  plasma  and  a third  luminous  region 
against  the  Teflon  surface.  At  5 u s the  discharge  current  has  just  reached  its  negative 
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Figure  5.  Simultaneous  High  Speed  Photographs  of  Side  & Front  of  Aeeelerator  - 
Side  View  on  Left,  Anode  Left  - Front  View  on  Right  - Anode  Bottom 
Time  in  Microseconds  at  Left  of  La  eh  Photograph. 
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peak.  The  second  luminous  plasma  continues  to  move  forward  while  the  third  appears 
to  remain  in  position  at  the  surface  of  the  propellant.  A close  scrutinty  of  the  photo- 
graph at  8 us  indicates  that  the  third  luminous  plasma  has  separated  from  the  ablator 
surface.  The  second  positive  maximum  discharge  current  occurs  at  roughly  7.8  us 
or  just  prior  to  the  third  arc's  separation.  Beyond  this  point  in  time  the  arc  becomes 
insufficiently  luminous  to  draw  conclusions  from  the  photographs.  The  images  are 
not  as  clearly  defined  and  are  barely  evident  when  printed. 

Certain  conclusions  concerning  the  evolution  of  accelerated  plasma  can  be  drawn 
from  the  above  experimental  evidence.  It  is  first  obvious  that  the  plasma  is 
accelerated  in  discrete  bunches  which  are  called  plasmoids.  The  appearance  of 
these  plasmoids  seems  to  be  coincident  with  each  quarter  cycle  of  the  current  and, 
furthermore,  each  plasmoid  appears  to  begin  accelerating  away  from  the  ablator 
surface  at  or  slightly  after  a maxima  of  the  discharge  current.  It  was  shown, 
however,  through  an  analytic  treatment  of  the  one  dimensional  arc  dynamics,  that 
the  separation"  of  each  plasmoid  off  the  ablation  surface  actually  becomes  incipient 
at  the  instant  the  direction  of  the  gradient  of  magnetic  flux  density  reverses.  This 
reversal  in  the  gradient  of  B at  the  ablator  surface  is  the  result  of  the  change  in 
direction  of  current  flow  at  the  ablator,  which  leads  the  reversal  of  total  discharge 
current  by  a small  amount  of  time  because  of  the  inductance  between  the  capacitor 
bank  terminals  and  the  point  on  the  electrodes  at  which  the  ablator  is  positioned. 

Thus,  each  plasmoid  begins  to  separate  from  the  ablator  surface  just  prior  to  each 
reversal  of  the  discharge  current. 

There  is  no  apparent  discontinuity  in  the  plasma  motion  at  the  current  reversal  times 
even  though  the  total  discharge  current  is  instantaneously  zero  at  these  times.  No 
so-called  "restrikes"  are  obvious  from  these  photographs.  This,  and  other  supporting 
data,  led  to  the  hypothesis  that  each  subsequent  plasmoid  is  in  reality  a portion  of  the 
previous  standing  arc  which  has  been  "sheared"  off  by  the  action  of  strong  magnetic 
pressure  forces  generated  once  the  direction  of  current  flow  at  the  ablator  surface 
has  changed.  The  continuous  depolymerization  of  Teflon  allows  a sufficient  quantity 
of  new  mass  to  enter  the  interelectrode  gap  to  keep  the  arc  sustained  through  the  zero 
total  discharge  current  points.  Moreover,  the  current  flow  through  each  subsequent 
standing  arc  is  in  part  also  flowing  through  the  accelerating  arc,  thereby  coupling  the 
two  while  they  coexist  within  the  interelectrode  gap. 
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With  regard  to  the  two  dimensional  shape  of  the  arc,  it  is  difficult  to  ascertain  the 
true  arc  configuration  from  these  photographs  since  the  luminosity  is  very  intense. 

The  arc  docs  appear  to  be  fairuly  uniform  across  the  gap  up  until  about  3 microseconds 
but  beyond  that  one  can  detect  a certain  degree  of  nonuniformity,  especially  obvious 
in  the  photographs  at  4.5,  5. 0,  and  5.5  microseconds.  The  difficulty  in  interpreting 
the  true  arc  shape  from  these  photographs  prompted  the  use  of  time  of  arrival  probing 
to  determine  the  leading  edge  shape  and  the  result  of  these  measurements  are  reported 
upon  next. 

3.3  Arc  Velocity  and  Shape  Using  Time-of- Arrival  Probe 

A double  electrostatic  probe  was  constructed  using  0.125  inch  (3.2  mm)  diameter 
semi-rigid  coaxial  cable.  The  outer  and  inner  conductors  were  the  cathode  and  anode 
of  the  probe;  respectively.  A 9V  battery  was  used  to  drive  the  probe  and  the  output 
signal  was  detected  using  a current  transformer.  The  outside  surface  of  the  outer 
conductor  was  insulated  with  a coating  of  high  temperature  epoxy. 

The  current  tansformer  output  and  total  discharge  current  were  displayed  simultaneously 
on  a dual  beam  analog  oscilloscope.  The  discharge  current  provided  a common  reference 
for  time  zero,  needed  to  do  the  data  reduction.  Several  oscillograms  were  recorded 
at  each  location  probed  and  these  data  were  averaged.  Due  to  significant  uncertainty, 
it  was  not  possible  to  decifer  accurate  time-of-arrival  any  closer  to  the  ablator  sur- 
face than  about  3. 5 cm.  The  maximum  distance  allowed  by  the  facility  from  the 
ablator  surface  to  the  probe  was  12  cm,  or  roughly  3 cm  beyond  the  exit  plane  of  the 
electrodes. 

In  order  to  determine  the  shape  of  the  front  at  given  times  the  raw  data  on  time-of- 
arrival  of  the  first  plasmoid  was  plotted  as  a function  of  axial  distance  at  each  lateral 
position  between  the  electrodes.  A curve  was  then  drawn  through  the  data  points  and 
the  position  of  the  front  at  the  desired  time  was  read  from  each  of  these  plots.  These 
points  were  then  plotted  as  a function  of  lateral  position  between  the  electrodes.  The 
resulting  curves  are  plotted  in  Figure  6 which  shows  the  shape  of  the  accelerating 
front  at  two  different  times.  The  line  x - o corresponds  to  the  ablator  surface  and 
the  line  y = o is  the  anode  surface. 
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The  first  conclusion  drawn  from  the  data  is  that  the  front  is  moving  faster  at  the 
anode  than  at  the  cathode.  Indeed,  the  raw  time-of-arrival  data  indicates  that  the 
front  velocity  at  the  cathode  is  approximately  one-half  that  at  the  anode.  Further- 
more, the  front  appears  to  be  fairly  uniform  over  between  50  and  75%  of  the  inter- 
electrode gap  from  the  anode  up.  The  raw  data  also  shows  that  the  uniform  portion 
of  the  front  is  moving  at  fairly  constant  velocity  after  accelerating  during  approxi- 
mately the  first  1.5  microseconds  and  that  the  velocity  near  the  cathode,  while 
increasing  up  to  about  2 microseconds,  begins  to  slowly  decrease  beyond  this  time. 

No  conclusive  information  concerning  the  second  or  other  subsequent  plasmoids 
could  be  drawn  from  the  raw  time-of-arrival  data.  While  a second  discrete  peak  in 
the  output  signal  was  discernible  in  some  of  the  oscillograms,  this  peak  appeared 
to  be  merged  with  the  first  in  some  cases  and  appeared  to  be  nonexistent  in  others. 

It  was  not  understood  why  this  was  the  case,  but  it  was  presumed  that  it  was  due  to 
the  statistical  nature  of  such  probes  in  detecting  the  presence  of  a plasma  via  the 
presence  of  its  conductivity.  This  explanation  seemed  plausible  since  the  shape  of 
the  peaks  even  at  a given  position  were  not  identical  from  shot-to-shot  but,  on  the 
contrary,  differed  quite  a bit  from  each  other,  even  though  the  onset  of  current  flow 
(time  of  arrival)  was  fairly  reproducible. 

One  of  the  more  interesting  observations  concerning  this  data  is  the  absence  of  any 
apparent  boundary  layer  at  the  electrode  surfaces.  Even  though  the  closest  point 
at  which  the  probe  could  be  positioned  to  the  electrodes  was  3 mm  off  the  surface, 
at  a distance  of  10  cm  downstream  from  the  ablator  surface  one  would  expect  to  see 
some  boundary  layer  structure  at  this  distance  above  the  electrode  surface.  One 
fact  is  certain  on  the  basis  of  front  location  and  that  is  that  the  presence  of  viscosity 
in  a continuum  sense  is  not  a significant  factor  in  controlling  the  motion  of  the  plasma 
front  near  the  electrode  surfaces.  This  is  probably  the  case,  at  least  with  regard  to 
the  motion  of  the  front,  because  the  plasma  density  at  the  front  is  low  and  the  no-slip 
condition  at  the  electrodes  is  probably  not  applicable  in  the  vicinity  of  the  leading 
edge. 

3.4  Distribution  of  Magnetic  Flux  Density 

Small  coil  magnetic  probes  were  constructed  for  measurement  of  magnetic  flux 
density  within  the  interelectrode  gap.  The  probes  were  wound  using  . 008  inch  (0. 2 mm) 
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diameter  lacquered  copper  wire.  Each  probe  consisted  of  9 complete  closely  spaced 
turns  on  an  I.  D.  of  .069  inch  (1.75  mm).  The  ends  of  the  probes  were  soldered  to 
miniature  coaxial  cable.  Dense  alumina  tubing  having  anO.  D.  of  6 mm  and  I.  D.  of 
3 mm  was  used  to  house  the  probe  and  coaxial  line.  The  end  of  the  tubing  containing 
the  probe  was  potted  to  a hemispherical  shape  with  high  temperature  resistant  epoxy. 
The  epoxy  potting  was  necessary  to  protect  the  probe  from  the  hot  plasma  and  was 
advantageous  since  it  allowed  a more  streamlined  shape  to  be  molded  at  the  tube  end. 
A photograph  of  a magnified  probe  winding  is  presented  in  Figure  7. 

Data  were  taken  on  the  Nicolet  Model  1090A  digital  memory  oscilloscope.  The  probe 
was  moved  in  .25  cm  increments  in  both  directions,  yielding  excellent  spatial  resolu- 
tion. The  probes  were  calibrated  using  a Helmholtz  coil  and  showed  flat  response 
iqj  to  10  MHz. 

In  our  first  year’s  work,  measurements  of  magnetic  flux  density  only  along  the 
accelerator  centerline  were  made.  It  was  seen  there  that  the  temporal  variation  of 
B at  the  ablator  surface  followed  the  discharge  current  to  within  a few  percent  over 
the  entire  cycle.  The  analytic  one-dimensional  prediction  of  the  appropriate  con- 
stant of  proportionality  between  B and  i indicated  an  arc  width  of  2. 93  cm  which  was 
in  surprisingly  good  agreement  with  the  luminous  arc  width  measured  off  of  the  high 
speed  photographs.  Following  up  with  this  in  two  dimensions,  the  ratios  i/imax  and 
B/Bniax  at  various  locations  along  the  surface  of  the  ablator  were  plotted  simultan- 
eously as  a function  of  time.  This  plot  is  presented  in  Figure  8.  The  first  observa- 
tion to  be  made  from  these  data  is  that  the  magnetic  flux  is  very  nearly  proportional 
to  the  discharge  current  everywhere  along  the  surface  of  the  ablator,  but  there  is 
some  distinction  between  the  first  and  second  half  cycles  of  the  current.  During  the 
first  half  cycle  there  is  a slight  gradient  in  B along  the  surface  of  the  ablator  such 
that  B is  increasing  from  roughly  the  center  of  the  interelectrode  gap  towards  the 
electrodes.  The  second  half  cycle  is  characterized  by  a gradient  in  B with  B 
decreasing  from  the  cathode  to  the  anode.  The  minimum  negative  value  of  B at  the 
anode  is  -0. 50T  and  at  the  cathode  it  is  -0. 81 T.  The  slight  gradient  evident  in  the 
first  half-cycle  is  within  the  realm  of  experimental  error,  but  there  is  no  question 
about  the  fact  that  the  observed  gradient  during  the  second  half  cycle  is  real. 

, 
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Figure  8.  Magnetic  Probe  Response  and  Discharge  Current 
versus  Time  at  Various  Positions  Along  the  Ablator  Surface 

A gradient  in  the  direction  indicated  would  imply  a finite  flow  of  current  away  from 
the  surface  of  the  ablator  (i.  e. , in  the  axial  direction)  as  indicated  by  Amperes 
law  in  differential  form  which  gives 

J1  = Bj/M  J2  = -V*  (2) 

where  and  J ^ are  the  current  density  components  in  the  x-  and  y-  directions; 
respectively,  and  subscripts  x and  y indicate  differentiation  in  those  directions. 

A curvature  in  the  lines  of  current  flow  away  from  the  ablator  is  dictated  by  the 
above  equations  and  the  observed  direction  of  By.  This  curvature  was  more  or  less 
expected  to  exist  once  the  time-of-arrival  data  had  been  obtained  since  those  data 
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showed  that  the  acceleration  front  was  moving  faster  at  the  anode  surface  than  at 
the  cathode  surface,  implying  a thicker  arc  at  the  anode  as  compared  to  tlu*  cathode. 

The  above  observations  are  more  clearly  indicated  in  Figure  9 where  the  magnetic 
flux  lines  have  been  plotted  at  four  stages  of  arc  progression  during  the  discharge. 

The  first  plot  (Figure  9a)  shows  the  lines  of  constant  B at  the  instant  of  peak  current 
in  the  first  half-cycle  (i.e. , approximately  1.93  u s from  initiation).  The  lines 
begin  at  a value  of  1.15T  closest  to  x = o and  the  increment  between  lines  is  -0.165T. 
At  this  instant  of  time  the  presence  of  a standing  arc  is  clearly  evident  from  the  very 
rapid  decrease  of  B from  the  value  at  x = o.  It  is  also  obvious  that  the  y-component 
of  the  current  density  (J2)  is  considerably  larger  at  the  cathode  than  at  the  anode. 
Numerical  differencing  indicates  a gradient  of  .53T/cm  near  the  ablator  which  trans- 
lates to  a current  density  of  42  MA/ m2  at  the  cathode.  Similar  calculations  reveal 
a current  density  of  25  MA/m2  at  the  anode  surface  close  to  the  ablator.  Further 
analysis  of  the  data  indicates  that  the  current  density  is  fairly  constant  at  the  surface 
of  both  the  anode  and  the  cathode  out  to  a distance  of  approximately  1 cm  at  the  cathode 
and  2. 8 cm  at  the  anode,  after  which  the  current  density  falls  off  rapidly  until  the 
B-  field  reaches  its  minimum  negative  value.  At  that  point  the  current  density  is 
zero  and  one  would  suspect  that  this  would  define  the  maximum  downstream  extremity 
of  the  arc.  It  is  seen  that  this  is  truly  the  case  once  the  data  obtained  using  time-of- 
arrival  probes  is  superimposed  on  the  flux  density  contours.  The  leading  edge  posi- 
tion shown  in  Figure  6 has  been  superimposed  as  a dotted  line  on  Figures  9a 
and  9b.  In  Figure  9b  the  flux  lines  one  microsecond  after  maximum  current 
(2.93  n s after  initiation)  are  plotted.  At  this  instant  of  time  the  high  speed  photo- 
graphs do  not  indicate  arc  separation  and  this  is  supported  by  the  magnetic  probe 
data.  The  flux  line  closest  to  x = o in  Figure  9b  represents  a field  of  . 99T  with  an 
increment  of  -.  165T  between  subsequent  lines.  At  this  instant  the  cathode  current 
density  has  decreased  to  4.5  MA/m2  near  the  ablator  and  the  anode  current  density 
has  also  decreased  in  value  to  .67  MA/m2.  Again,  the  current  density  at  the  anode 
and  cathode  surfaces  appears  to  be  fairly  constant  up  until  a downstream  distance 
just  short  of  the  plasma  front.  Skewing  of  the  current  density  lines  is  more  obvious 
in  Figure  9b,  however,  as  can  be  interpreted  from  gradients  of  B in  the  y-direction. 

Perhaps  the  most  Interesting  flux  plot  occurs  at  the  instant  of  zero  current  (Figure  9c) 
which  occurs  3. 85  n s after  initiation.  At  that  time  the  magnetic  flux  still  shows  small 
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positive  values  at  the  ablator  surface  which  range  between  . 07T  and  . 2T,  but  the 
axial  gradient  of  B has  changed  sign.  The  lino  closest  to  the  ablator  surface  (x  o) 
has  a value  of  . 08;! T and  the  increment  between  lines  is  t.  083 T up  |o  the  central 
closed  loop  representing  .413T,  and  then  083T  beyond  that  closed  loop.  The 
observed  change  in  Bx  from  positive  to  negative  at  the  ablator  surface  implies  that 
the  current  density  at  the  wall  has  changed  direction.  This  is  of  interest  because 
the  total  discharge  current  has  not  yet  changed  direction  and,  in  fact,  is  zero  at  the 
instant  of  time  depicted  in  Figure  9c.  No  conclusions  can  be  made  regarding  the 
flow  of  current  in  the  plasma  at  this  instant  solely  on  the  basis  of  the  flux  density 
distribution,  but  the  measurements  of  local  current  flow  using  a current  probe  has 
shed  considerable  light  on  the  situation  and  those  results  appear  in  the  next  section. 

It  is  obvious  from  Figure  9c,  however,  that  separation  of  the  first  plasmoid  has 
occurred.  This  is  clear  from  the  magnetic  flux  distribution  since  the  absolute  value 
of  the  axial  gradient  of  B has  two  local  maxima  within  the  plotted  demain.  These 
two  maxima  are  located  in  the  intervals  0 < x < 1 cm  and  6 cm  < x < 8. 5 cm.  The 
appearance  of  these  maxima  in  | Bj,  | manifest  the  existance  of  two  separate  arcs ; 
one  at  the  ablator  surface,  and  one  downstream  of  it.  Moreover,  the  relatively  large 
value  of  By  near  the  cathode  surface  in  particular  indicates  considerable  curvature 
in  the  lines  of  current  flow  near  the  cathode.  More  will  be  said  about  this  in  the 
next  section. 

Figure  9d  is  a plot  of  the  flux  lines  one  microsecond  after  zero  current  (4. 85  4 s 
after  initiation).  The  line  closest  to  the  ablator  has  the  value  -.495T  and  the  incre- 
ment between  lines  is  . 083T.  The  current  density  at  the  ablator  surface  near  the 
cathode  has  a value  of  -28.4  MA/m^  and  at  the  anode  surface  the  value  is  -10  MA/m2. 
Thus,  the  existence  of  higher  local  current  density  at  the  cathode  as  compared  to  the 
anode,  noted  during  the  first  half  cycle,  is  also  the  case  during  the  second  half  cycle. 
The  large  axial  gradient  in  B which  extends  from  the  ablator  surface  out  to  about 
1. 8 cm  defines  the  extent  of  the  second  plasmoid  at  the  instant  of  time  shown.  A 
second  local  large  gradient  of  B downstream  from  x = o is  no  longer  present  at  this 
instant  of  time,  which  leads  one  to  conclude  that  the  first  plasmoid  has  already  moved 
out  of  the  interelectrode  gap.  To  what  extent  this  observation  is  true  will  be  further 
analyzed  in  the  next  section  of  this  report. 
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No  further  analysis  of  the  data  will  be  made  beyond  4.  85  us  of  the  discharge  cycle, 
the  reason  for  this  being  that  a reasonably  accurate  interpretation  of  the  data  becomes 
difficult  after  this  time  due  to  the  relatively  low  level  of  the  measurements.  It  is 
noted,  however,  that  the  most  interesting  aspects  of  plasma  acceleration  in  this 
device  have  already  occurred  by  this  time.  Namely,  the  initial  generation  of  plasma; 
the  separation  of  the  standing  arc  from  the  ablator  surface  and  its  subsequent 
acceleration  out  of  the  interelectrode  gap;  and  the  regeneration  of  a second  standing 
arc  which  coexists  with  the  accelerating  arc  in  a rather  unique  way.  An  understanding 
of  these  phenomena  will  lead  to  considerable  insight  into  the  parameters  which  make 
ablative  plasma  acceleration  so  unique. 

3. 5 Current  Distribution 

In  order  to  obtain  some  idea  of  the  current  distribution  within  the  plasma  a Rogowski 
loop  was  constructed  which  could  be  positioned  at  known  locations  from  the  propellant 
surface  along  the  centerline.  Since  the  lc^p  would  be  immersed  in  the  hot  plasma  it 
was  necessary  to  protect  the  winding  by  potting  the  entire  coil  with  a high  temperature 
epoxy  after  it  was  wound.  The  loop  was  rectangular  in  shape,  having  an  inside  width 
and  length  of  2 mm  and  10.7  cm.  The  coil  was  designed  so  as  to  extend  laterally 
slightly  beyond  the  ends  of  the  electrodes  to  insure  that  all  current  flowing  through 
the  plasma  at  a given  location  would  be  intercepted.  The  diameter  of  the  winding, 
including  the  opoxy  coating,  was  2.8  mm.  The  coil  was  calibrated  on  the  same  circuit 
as  that  used  to  measure  the  total  discharge  current. 

Data  were  obtained  at  each  position  scanned  by  discharging  the  capacitor  bank  and 
recording  the  total  discharge  current  and  the  plasma  current  probe  response  simultane- 
ously on  the  digital  memory  oscilloscope.  The  total  discharge  current  was  used  to 
insure  a consistent  zero  time  reference  for  all  data.  It  is  noted  that  with  the  coil 
immediately  up  against  the  Teflon  surface  the  center  of  the  loop  was  3.6  mm  away 
from  it.  The  current  flowing  at  the  ablation  surface  was  obtained  by  extrapolating 
back  using  the  two  data  points  closest  to  it. 

Figure  10  shows  the  distribution  of  current  at  four  different  times  during  the  discharge. 
These  four  instants  of  time  are  the  same  as  those  at  which  magnetic  flux  was  plotted. 
Figure  10a  shows  the  distribution  of  current  flowing  through  the  plasma  at  1.93  u s 
which  corresponds  to  the  instant  at  which  peak  discharge  current  occurs.  The  data  is 
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presented  in  the  form  of  percentage  contours,  the  line  closest  to  x = o representing 
the  farthest  downstream  position  at  which  20%  of  the  total  current  flowing  at  that  time 
is  flowing  through  the  plasma,  and  each  subsequent  line  showing  10%  increments  of 
the  total.  The  same  presentation  appears  in  Figure  10b  at  2.93  y s. 

It  is  clear  from  Figure  10a  that  the  discharge  current  is  distributed  fairly  evenly 
throughout  the  arc  in  each  plane  parallel  to  the  electrodes  although  significant  variations 
in  arc  thickness  are  evident  from  one  plane  to  the  next.  The  arc  is  thicker  at  the 
cathode  than  at  the  anode,  which  is  what  estimates  of  the  current  density  from  the 
magnetic  flux  seemed  to  indicate.  Having  these  two  sets  of  data  it  is  possible  to 
calculate  an  average  width  of  the  arc  at  both  the  anode  and  the  cathode.  Calculation 
of  J2  using  Equation  (2)  and  the  magnetic  flux  data  yieleded  values  of  42  MA/m^  and 
25  MA/m2  at  the  cathode  and  anode;  respectively  near  the  ablation  surface.  The 
current  distribution  indicates  that  near  the  cathode  20%  of  the  29.9  kA  is  flowing  out 
to  approximately  0.  2 cm  and  near  the  anode  this  same  amount  of  current  is  flowing 
out  to  0.7  cm.  These  values  lead  to  an  estimate  of  0. 7 cm  for  the  arc  width  at  the 
cathode  and  3.4  cm  at  the  anode. 

At  2. 93  jis  (Figure  10b)  the  current  is  still  fairly  uniformly  distributed  throughout  the 
arc.  Estimates  of  the  arc  width  based  on  data  at  this  time  level  are  approximately  the 
same  as  those  values  generated  at  1.93  us.  The  degree  to  which  the  arc  leading  edge 
at  the  anode  is  preceding  the  leading  edge  at  the  cathode  is  considerably  more  severe 
at  this  level  of  time,  however.  It  is  noted  that  the  arc  at  the  cathode  has  not  advanced 
downstream  at  all  during  the  one  microsecond  between  Figures  10a  and  10b,  whereas 
there  has  been  a movement  of  approximately  2.5  cm  at  the  anode  surface  during  the 
same  time  period.  Comparison  of  the  data  in  Figures  l0.a  and  10  b with  the  plasma 
front  location  determined  using  the  time-of-arrival  probe  indicates  that  at  1.93  u s 
the  plasma  front  and  arc  leading  edge  are  very  close  to  each  other  at  the  anode  and 
cathode  surfaces  but  within  the  gap  the  plasma  leading  edge  is  ahead  of  the  arc  front 
as  much  as  1 cm  at  some  locations.  At  2.93  u s the  situation  is  somewhat  different, 
however.  The  plasma  leading  edge  at  the  cathode  surface  is  approximately  2.2  cm 
ahead  of  the  arc  front  and  at  the  anode  the  plasma  leading  edge  precedes  the  arc  front 
by  some  1.3  cm.  Within  the  gap  one  finds  that  the  plasma  leading  edge  is  a maximum 
of  2. 8 cm  ahead  of  the  arc  front  at  y = 5 cm.  These  observations  indicate  that  the 
plasma  is  indeed  accelerating  out  ahead  of  the  arc  even  as  the  arc  itself  spreads 
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downstream.  This  is  attributed  to  the  action  of  the  Lorentz  force  on  charged 
particles  in  the  arc  front.  It  is  also  noted  that  the  locus  of  points  at  which  the 
magnetic  field  is  a minimum  negative  value  corresponds  quite  well  to  the  arc  front 
location  as  would  be  suspected.  Hence,  those  charged  particles  which  are  accelerated 
out  ahead  of  the  arc  will  find  themselves  in  a favorable  crossed  electromagnetic  field 
and  will,  therefore,  continue  to  accelerate. 

The  data  at  3.  85  u s (discharge  current  reversal)  is  presented  in  a slightly  different 
way  because  no  current  is  flowing  from  the  capacitors  at  that  instant.  Current  is, 
however,  flowing  within  the  plasma  but  conservation  of  charge  requires  that  no  net 
current  flow  exist.  This  is  indeed  the  case  since  the  reversal  of  current  flow  direction 
within  the  plasma  occurs  roughly  one  half  microsecond  prior  to  total  discharge  current 
reversal  and  begins  flowing  in  the  opposite  direction  at  the  ablation  surface.  In 
Figure  10c,  therefore,  the  dashed  line  shows  the  maximum  downstream  position  up  to 
which  the  current  is  reversed.  The  solid  contours  on  either  side  indicate  the  extremities 
of  various  percentages  of  the  total  reverse  current  flow.  Between  x = o and  the  dashed 
line,  the  first  solid  contour  shows  50%  of  the  total  and  each  subsequent  line  includes 
10%  more.  To  the  right  of  the  dashed  line  the  first  solid  contour  shows  50%  and  each 
subsequent  line  10%  less  net  current  flow.  The  two  separate  distinct  arcs  which  are 
formed  as  a result  of  the  current  reversing  its  direction  of  flow  within  the  plasma  prior 
to  the  reversal  of  total  discharge  current  are  clearly  obvious  in  Figure  10c.  Moreover, 
because  of  the  close  proximity  of  the  two  arcs  near  the  cathod  surface,  it  appears  that 
there  is  indeed  a closed  current  loop  within  the  plasma  at  the  time  no  net  discharge 
current  is  flowing  from  the  capacitor  bank.  The  arc  at  the  ablation  surface  is  being 
sustained  at  this  instant  of  time  by  energy  previously  deposited  in  the  accelerating  arc. 

In  Figure  lOd,  one  microsecond  after  discharge  current  reversal,  the  data  are  again 
presented  as  percentage  of  total  discharge  current  at  that  instant,  the  first  contour 
representing  20%  and  each  subsequent  contour  10%  more.  This  figure  indicates  that 
the  accelerating  arc  has  moved  entirely  out  of  the  electrode  gap  and  shows  the 
existence  of  the  developed  second  standing  arc. 

A point  worth  clarification  with  regard  to  the  separating  or  degeneration  of  the  initial 
arc  into  two  separate  arcs  is  the  direction  of  accelerating  forces  just  prior  to  this 
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occurrence.  Let  us  trace  the  events  from  the  time  peak  current  occurs  and  consider 
only  the  Lorentz  force  on  the  plasma,  which  is  proportional  to  the  x-derivative  of  the 
* square  of  the  magnetic  flux  density  (when  9 B2/9  x is  negative  the  Lorentz  force  is 

positive).  Maximum  flux  density  at  the  wail  occurs  at  maximum  discharge  current, 
as  observed  from  the  data  on  B.  Moreover,  the  values  of  B at  the  ablator  are  the 
local  maximum  values  at  any  distance  y above  the  anode.  Hence,  the  gradient  in 
magnetic  pressure  is  negative  up  to  the  time  peak  current  occurs  and,  therefore, 
the  Lorentz  force  is  directed  away  from  the  ablator,  tending  to  accelerate  the  plasma. 
Consider  now  what  happens  as  the  discharge  current  decays  from  the  maximum. 

The  magnetic  flux  density  at  the  ablation  surface  will  also  decrease  from  its  maximum 
value.  Due  to  the  conductivity  of  the  plasma,  however,  particles  which  have  been 
ablated  up  until  that  time  are  carrying  magnetic  field  energy  with  them  away  from  the 
i wall,  and  as  these  particles  move  further  downstream  the  conversion  of  magnetic 

field  energy  into  kinetic  energy  via  the  Lorentz  force  will  cause  the  amount  of  magnetic 
field  energy  associated  with  each  accelerating  particle  to  decrease.  The  rate  of 
decrease  of  magnetic  field  energy  within  a small  region  close  to  the  ablation  surface 

t 

( and  the  rate  of  decrease  of  magnetic  flux  density  at  the  surface  itself  arc  the  relevant 

parameters  in  determining  what  happens  next.  If  these  rates  are  comparable  then  the 
gradient  in  B2  at  the  ablation  surface  will  be  small  but  favorable  regarding  accelera- 
tion, but  as  soon  as  the  rate  of  decrease  in  B at  the  ablation  surface  is  faster  than 

I 

the  rate  of  energy  conversion  from  magnetic  to  kinetic  within  the  small  volume  close 
to  the  wall,  the  gradient  of  B will  become  adverse.  Simultaneously  the  electric  field 
at  the  ablation  surface  must  change  direction  and  the  direction  of  current  flow  at  this 
surface  must  reverse.  The  most  interesting  aspect  of  this  is  the  fact  that  the  total 
discharge  current  flowing  from  the  capacitor  has  not  changed  its  direction  of  flow 
at  the  time  the  gradient  in  B2  has  become  adverse  and  hence,  the  vast  majority  of 
the  plasma  produced  up  until  that  time  is  still  being  accelerated  with  a favorable 
gradient  in  B2  present  within.  The  arc  is,  therefore,  forced  to  split  into  two  separate 
arcs,  one  carrying  current  in  the  same  direction  as  the  current  flow  from  the  capacitor, 
and  one  carrying  current  in  the  opposite  direction.  This  leads  to  a closed  current  loop 
within  the  plasma  which  amplifies  the  magnetic  field  strength  between  them  and  pro- 
duces as  a result  an  "explosion"  force  distribution,  tending  to  push  the  accelerating  arc 
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out  of  the  electrode  gap  and  also  tending  to  push  the  newly  formed  reverse  current 
arc  against  the  ablation  surface.  This  can  be  more  clearly  visualized  once  it  is 
realized  that  the  direction  of  current  flow  at  the  ablation  surface  will  reverse  while 
the  magnetic  flux  density  is  still  positive  and  this  implies  a negative  Lorenlz  force 
on  the  particles  carrying  the  current.  This  situation  persists  only  up  until  the  time 
the  total  discharge  current  reverses,  at  which  time  the  magnetic  flux  density  also 
changes  direction,  and  the  forces  on  the  standing  arc  become  favorable  regarding 
its  acceleration.  As  nearly  as  can  be  estimated,  the  condition  of  negative  Lorenlz 
force  on  the  second  standing  arc  exists  for  approximately  one  half  microsecond 
prior  to  discharge  current  reversal.  That  is,  arc  degeneration  becomes  incipient 
at  about  3.3  us  subsequent  to  initiation. 

3, 6 Electron  Number  Density  and  Temperature 

3.6. 1 Description  of  Langmuir  Probe 

In  order  to  measure  electron  number  density  and  temperature  within  the  accelerating 
plasma  a conventional  double  Langmuir  probe  was  constructed.  The  probe  consisted 
of  two  strands  of  5 mil  (.  126mm)  tungsten  wire  enclosed  within  a glass  rod.  The 
length  of  wire  exposed  beyond  the  end  of  the  glass  rod  was  1 cm  (see  Figure  11). 


Figure  11.  Langmuir  Probe  Configuration 

Initial  experimentation  with  this  probe  led  to  an  unexpected  problem  associated  with 
its  utilization  in  the  plasma  produced  from  depolymerized  Teflon.  It  was  found  that 
a tendency  for  the  tungsten  wires  to  become  contaminated  by  the  plasma  was  prevalent. 
Moreover,  this  contamination  was  severe  enough  to  make  acquisition  of  a single 
voltage/current  characteristic  impossible.  A technique  for  remote  cleaning  of  the 
probe  wires  would  therefore  have  to  be  developed  if  experimentation  was  to  continue. 
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The  technique  finally  decided  upon  was  to  use  electron  bombardment  of  the  probe 
wires.  The  electron  source  for  this  bombardment  was  a tungsten  ribbon  5 mils 
(.  127mm)  thick  and  20  mils  (.508mm)  wide  wound  in  a helical  fashion  with  three* 
complete  turns  on  a diameter  of  roughly  5 cm.  The  probe  was  |x>sitioned  at  the 
geometrical  center  of  this  winding  and  the  ends  of  the  ribbon  were  fastened  to  two 
standoffs  fixed  to  the  body  of  the  probe.  Experimentation  with  this  arrangement 
resulted  in  the  appropriate  parameters  necessary  to  do  a good  cleaning  job  on  the 
probe  wires.  It  was  found  that  by  running  approximately  230W  at  33  VAC  through  a 
transformer  directly  to  the  tungsten  ribbon  and  at  the  same  time  maintaining  a DC 
voltage  of  1900  to  2300V  on  the  probe  wires,  a current  of  about  70MA  would  pass 
from  the  ribbon  to  the  probe,  which  was  sufficient  to  heat  the  probe  wires  to  red 
hot  in  10  seconds.  The  final  steady  state  temperature  of  the  tungsten  ribbon  was 
2400°C,  measured  using  an  optical  pyrometer.  The  probe  voltage,  once  fully 
heated,  dropped  to  1000  VDC.  The  circuit  used  to  clean  the  probe  wires  is  depicted 
in  Figure  12. 


. 


33  VAC 


TUNGSTEN  RIBBON 
EMITTER 


70  mA 


LANGMUIR  PROBE 
1900  - 2300  VDC 


Figure  12.  Schematic  of  Probe  Cleaning  Circuitry 
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A second  a 17  result  of  cleaning  the  probe  wires  was  local  destruction  of  the  glass  rod 
enclosure  near  the  point  from  which  the  wires  protruded.  In  order  to  preserve  the 
lifetime  of  a given  probe  this  localized  erosion  had  to  be  minimized  so  a second  glass 
enclosure  in  the  form  of  a 2cm  ID  tube  was  placed  around  the  glass  rod  such  that  the 
end  of  the  tube  was  just  ahead  of  the  end  of  the  glass  rod.  This  solved  the  problem 
of  local  destruction  at  the  end  of  the  glass  rod.  A photograph  of  the  final  Langmuir 
probe  configuration  with  remote  cleaning  capability  built  in  is  presented  in  Figure  13. 
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Figure  13.  Langmuir  Probe  Equipped  for  Remote  Cleaning 
by  Electron  Bombardment 


3.6.2  Data  Acquisition  Procedure 

In  order  to  measure  electron  number  density  and  temperature  at  a given  point 
within  a plasma  it  is  necessary  to  obtain  a complete  voltage/current  characteristic 
at  that  point.  A good  characteristic  has  the  shape  of  an  arc  tangent  function, 
asymptoting  the  ion  saturation  values  of  current  as  the  applied  voltage  is  made 
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large  positive  or  negative,  and  crossing  through  the  origin  of  the  IV-plane.  Under 
conditions  of  a plasma  which  is  either  quiescent  or  moving  at  a constant  speed  one 
can  normally  pulse  the  probe  with  a voltage  ramp  function  of  short  duration  and 
display  I vs.  V at  the  probe  on  an  oscilloscope  having  appropriate  frequency  response 
In  a transient  situation,  however,  use  of  the  pulsed  probe  technique  would  require 
considerable  developmental  effort  to  manufacture  a driver  to  generate  a voltage  ramp 
of  sufficient  magnitude  in  the  short  time  period  necessary  to  obtain  a reasonable 
degree  of  resolution  (i.e.;  typically  -20V  to  +20V  in  approximately  1 microsecond). 
The  development  of  such  a driver  was  beyond  the  scope  of  the  current  program  so 
an  alternate  method  of  data  acquisition  was  formulated. 

The  technique  used  on  this  program  consisted  of  locating  the  probe  at  the  desired 
position,  applying  a DC  voltage  of  known  magnitude  across  it,  and  firing  the  accel- 
erator. Current  drawn  from  the  DC  source  as  the  plasma  passed  the  probe  was 
detected  using  a precalibrated  Rogowski  coil.  Discharge  current  and  the  probe 
current  were  simultaneously  displayed  on  a Tektronix  dual  beam  oscilloscope.  The 
discharge  current  provided  a common  time  base  to  which  all  data  could  be  referenced 
A new  oscillograph  is  recorded  for  each  DC  voltage  applied  across  the  probe  and 
reduction  of  the  data  to  a voltage/current  characteristic  at  a given  time  subsequent 
to  arc  initiation  is  accomplished  by  converting  the  probe  signal  trace  to  current  at 
that  instant  of  time  for  each  of  the  applied  voltages.  Typically,  up  to  ten  different 
DC  voltages  are  applied  at  any  given  location  in  order  to  generate  the  characteristic. 
A typical  volt-ampere  characteristic  obtained  as  outlined  above  is  presented  in 
Figure  14. 

3.6.3  Data  Reduction 

The  electron  number  density  and  temperature  are  related  to  the  volt-ampere 
characteristic  by  simple  probe  theory  through  the  following  two  relationships: 
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W \ "eAe(ll,606kTe/M+)5 

where  I is  current,  V is  voltage,  Iga(.  is  the  ion  saturation  current  (assumed 
symmetrical),  Tg  is  the  electron  temperature  in  ev,  ng  is  the  electron  member 
density,  A is  the  probe  collector  area,  e is  the  electron  change,  k is  the  Boltzman 
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constant  and  M is  the  ion  mass.  In  this  particular  case  the  collector  area  is 

o + 

.0157  cm  and  the  ion  mass  is  estimated  as  being  the  equivalent  of  a 2C  + 4F  gas 
mixture  or  16. 67amu.  Thus; 

n = 1.652  x 1015  I ,/T  2 
e sat  e 

with  n expressed  per  cubic  centimeter,  T in  ev  and  I in  amperes. 

6 0 Sell 

3.6.4  Number  Density  and  Temperature  Distribution 

The  modifications  made  to  the  probe  as  outlined  in  section  3.6. 1 in  order  to  keep 
the  probe  collectors  clean  made  it  impossible  to  scan  across  the  interelectrode 
gap  between  the  anode  and  cathode  because  of  space  limitations.  It  was,  however, 
possible  to  scan  along  the  centerline  in  the  axial  direction  but  the  degree  to  which 
this  could  be  accomplished  was  curtailed  by  the  fact  that  we  could  not  get  closer  to 
1 the  ablator  surface  than  5 cm.  The  reason  for  this  was  the  severe  electrical  noise 

generated  in  the  probe  circuitry  by  the  arc  initiation.  The  results  of  measurements 
made  along  the  accelerator  centerline  are  presented  in  Figure  15  where  electron 
number  density  as  a function  of  distance  is  plotted  at  various  times.  These  curves 
indicate  a continuous  monotonic  decrease  in  electron  number  density  with  down- 
stream distance.  The  overall  level  of  the  number  density  is  increasing  with  time 

as  the  arc  progi’esses  downstream.  The  maximum  value  of  the  number  density  which 

15  3 

was  measured  is  1.46  x 10  /cm  . If  one  were  to  linearly  extrapolate  this  value 

15  3 

back  to  the  surface  of  the  ablator,  a number  density  of  2.62  x 10  /cm  would  be 

obtained. 

1 i 

The  presence  of  conducting  plasma  as  a function  of  time  at  a given  point  can  be 

M 

more  or  less  interpreted  from  the  Langmuir  probe  response  at  that  point.  A 

atypical  oscillograph  of  this  response,  as  well  as  the  discharge  current  has  been 

reproduced  in  Figure  16.  Note  the  appearance  of  a very  clearly  defined  region  of 
conduction  between  6 and  10  microseconds  after  discharge  initiation.  The  probe 
response  during  this  interval  is  more  than  likely  due  to  the  passage  of  the  first 
plasmoid.  One  can  also  Interpret  the  peak  at  14  microseconds  as  being  due  to  the 
arrival  of  the  second  plasmoid. 

The  electron  temperature  along  the  centerline  of  the  accelerator  was  found  to  be 
constant  within  the  accuracy  of  the  experimental  results.  A value  of  lOev  is  character- 
istic and  did  not  vary  in  space  or  time.  This  result  is  not  surprising  but  is  mean- 
ingful because  it  indicates  that  the  energy  equipartition  time  for  electron/ion 
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Electron  Number  Density  as  a Function  of  Distance  Along  the 
Accelerator  Centerline  at  Three  Different  Times 
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Figure  16.  Typical  Langmuir  Probe  Response  and  Discharge 
Current  2cm  Past  the  Electrode  Exit  Plane 
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thermalization  is  much  longer  than  the  characteristic  flow  time  in  this 
accelerator. 

A two  dimensional  scan  across  the  interelectrode  gap  was  also  performed  at  a 
distance  of  2cm  beyond  the  exit  plane  of  the  electrode.  Space  limitations  within  the 
test  facility  limited  the  region  scanned  to  approximately  2cm  on  either  side  of  the 
centerline.  The  contours  as  a function  of  lateral  position  at  three  different  times 
are  presented  in  Figure  17.  These  contours  all  show  a local  maximum  in  electron 
number  density  near  the  centerline  of  the  accelerator  and  a rapidly  decreasing  level 
as  the  cathode  is  approached.  As  the  anode  is  approached  the  electron  number 
density  first  decreases  and  is  then  seen  to  be  increasing.  This  result  is  as  one 
would  normally  suspect  in  an  arc  of  this  type  except  for  the  local  minimum  at  5cm 
which  is  probably  more  related  to  the  flow  dynamics  than  to  the  distribution  of 
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Figure  17.  Electron  Number  Density  Across  the  Interelectrode 
Gap  a Distance  2cm  Past  the  Electrode  Fxit  Plane 
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4. 1 Physical  Model 

The  basis  for  the  physical  model  used  to  obtain  an  analytic  representation  of  the 
plasma  acceleration  mechanism  in  pulsed  ablative  thrusters  is  the  assumption  that 
the  plasma  can  be  considered  a continuum  at  all  points  in  space  and  for  all  time. 

Such  a model  results  in  a complete  description  of  the  plasma  by  specification  of  its 
macroscopic  thermodynamic  properties.  Changes  in  the  values  of  these  properties 
are  governed  by  a set  of  conservation  equations  which  are  derived  for  mass,  momen- 
tum, and  energy.  These  equations  are  valid  in  their  complete  form  for  the  solution 
of  any  flow  problem  and  only  the  boundary  and  initial  conditions  determine  precisely 
which  problem  is  being  solved.  The  consistency  and  accuracy  of  the  initial  and 
especially  the  boundary  conditions  is  therefore  essential  to  obtaining  an  analytic 
solution  to  the  governing  equations  which  accurately  describes  the  physical  process 
being  modeled. 

The  thermodynamic  state  of  the  plasma  is  assumed  to  be  a combination  of  ions, 
neutrals,  and  electrons,  each  of  which  is  in  local  thermodynamic  equilibrium  but 
thermodynamic  equilibrium  between  ions  and  electrons  is  not  assumed  to  exist  while 
ions  and  neutrals  are  assumed  to  be  in  equilibrium  with  one  another.  The  three 
species  are  considered  to  be  moving  at  the  same  baricentric  velocity  and  local 
charge  neutrality  is  assumed. 

4.1.1  Initial  Conditions 

A detailed  model  of  the  breakdown  or  "initiation"  phase  of  the  discharge  would  be 
required  if  one  desired  to  begin  calculation  of  the  plasma  flow  precisely  at  time 
t = 0.  Since  no  such  model  exists  at  this  time  for  the  breakdown  mechanism  used 
in  pulsed  ablative  thrusters  (i.e. , the  ignitor  plug  discharge  of  a microplasma  into 
the  interelectrode  gap)  it  is  not  possible  to  begin  calculation  at  precisely  t = 0. 
Instead,  calculations  begin  at  time  subsequent  to  t = 0 when  it  is  believed  that  the 
solid  Teflon  propellant  first  begins  to  depolymerize  into  gaseous  molecules.  This 
occurs  approximately  100  nanoseconds  after  Initiation,  ...  which  time  roughly  600A 
are  flowing  through  the  fully  developed  arc.  Up  until  this  time  the  arc  is  probably 
supported  mostly  by  the  microplasma  injected  by  the  ignitor  plug  and  to  some  degree 
by  consumed  electrode  material.  The  arc  thickness  at  the  time  calculations  are 
commenced  is  obtained  by  linear  extrapolation  of  the  early  experimental  data  to  that 
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time.  The  resulting  arc  thickness  obtained  in  that  fashion  is  2mm  and  the  arc  width  is 
3cm.  The  current  density  of  10  MA/m^  is  assumed  to  be  uniformly  distributed  from 
the  ablation  surface  to  the  downstream  extremity  of  the  arc.  The  magnetic  flux 
density  has  parabolic  decay  from  its  value  at  x = o to  its  value  at  the  leading  edge. 
These  values  are  computed  as  described  in  the  next  section.  The  initial  energy  per 
unit  volume  is  taken  as  the  energy  deposited  into  the  arc  up  until  100ns  divided  by  the 
volume  (4.57cm  ).  The  initial  arc  is  assumed  to  have  zero  velocity  and  uniform 
thermodynamic  pressure  and  mass  density,  obtained  by  subtracting  the  magnetic 
field  energy  from  the  initial  energy  input  and  assuming  a temperature  of  lOev. 
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4.1.2  Boundary  Conditions 

The  coordinate  system  is  located  with  its  origin  at  the  center  of  the  interelcctrode 
gap  on  the  ablator  surface  at  the  midpoint  of  the  ablator.  The  x-coordinatc  is  positive 
along  the  centerline  of  the  electrodes  with  positive  y directed  along  the  ablator  sur- 
face (toward  the  cathode)  and  positive  z parallel  to  the  electrodes.  The  assumption 
of  two  dimensional  motion  negates  variations  of  all  dependent  variables  in  the  z 
direction. 

The  mathematical  problem  is  to  find  the  solution  to  Equation  (3)  with  initial  conditions 
specified  in  the  domain  o^x^«  for  all  -^h  s y ^ + |h. 

The  existence  of  diffusion  terms  in  Equation  (3)  makes  them  a parabolic  set,  having 
only  one  characteristic  direction,  namely  streamlines,  given  by  dx/dt  = u and 
dy/dt  = v.  Since  the  direction  of  propogation  of  this  characteristic  will  always  be 
away  from  the  ablator  surface  (i.  e. , u ^ o always)  it  is  necessary  to  specify  a 
number  of  relationships  among  the  dependent  variables  equal  to  the  number  of 
dependent  variables  at  x = o for  all  y between  -h/2  and  -th/2.  Moreover,  since 
streamlines  will  always  propogate  away  from  the  electrode  surfaces  due  to  the 
probable  existence  of  boundary  layers,  all  dependent  variables  having  second  derivi- 
tives  in  the  y-direction  must  be  specified  on  y = ± h/2  for  all  x £ l . These  include 
B,  u,  v and  T.  Similarly,  these  dependent  variables  must  also  be  specified  ahead  of 
the  plasma  for  all  time. 

The  velocity  of  ablated  Teflon  at  the  instant  depolymerization  occurs  must  be  zero 
since  prior  to  that  time  the  gaseous  particles  will  not  have  experienced  any  imbalance 
of  forces  which  could  cause  them  to  accelerate.  Hence,  the  conditions  u = v = o 
must  exist  at  the  ablation  interface  for  all  time. 

The  temperature  at  the  ablation  surface  has  been  measured  under  controlled  thermal 
degradation  experiments.  It  was  found  that  the  temperature  of  the  solid  just 
prior  to  its  sublimation  is  not  a function  of  the  heat  flux  impingement  upon  it  and  is 
equal  to  640° C.  Thus,  a constant  temperature  surface  is  assumed  at  the  ablator. 

In  order  to  generate  a boundary  condition  on  the  density  it  is  necessary  to  resort  to 
a previously  published  empirical  data  correlation  between  the  total  mass  ablated 
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per  discharge  and  the  integral  of  the  square  of  the  discharge  current  over  time  16. 
It  was  shown  that  to  a fair  degree  of  accuracy  a significant  amount  of  data  obey  the 
relationship 


The  constant  a was  determined  precisely  for  the  experimental  apparatus  used 
throughout  this  program.  With  the  integral  of  current  squared  being  1.92  x 103A2s 
and  a mass  per  discharge  (averaged  over  1000  discharges  to  within  ± . 2u  g)  of 
0.148  mg,  the  constanta  is  7.7  x 10-3  g/A^s.  Extending  this  correlation  to  imply 
that  at  any  instant  of  time  t 

m(Ar)=  oc  ^ fdk 

leads  to  the  conclusion 

• * 2 

m = ocJL 

Furthermore,  by  assuming  a uniform  ablation  of  the  propellant  surface,  the 
instantaneous  mass  flow  rate  per  unit  area  can  be  obtained  by  dividing  m by  the 
propellant  area  Ap.  The  density  at  each  grid  point  on  the  ablation  surface  is 
calculated  using  a variation  of  the  finite  diffe  rence  scheme  outlined  in  the  Appendix. 
Letting  the  set  (i,  j)  represent  a point  on  the  ablator,  a control  volume  of  length 
Ax  and  height  2Ay  is  chosen  such  that  the  four  comers  of  this  volume  are  (i,  j + 1), 
(i*  j ~1)»  (i  + 1,  j + 1).  and  (i  +1,  j -1).  The  density  at  the  point  (i  + £,  j ) is  then 
calculated  to  second  order  accuracy  in  time  and  space  by  invoking  the  conservation 
of  mass  in  this  fixed  volume.  The  resulting  expressions  are 
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The  above  can  be  used  at  all  points  on  the  ablation  surface  except  those  points  in  the 
corners.  To  obtain  the  density  at  the  two  comer  points  values  at  the  two  adjacent 
points  on  the  ablation  surace  are  extrapolated. 

The  thermodynamic  pressure  at  the  ablation  surface  is  obtained  from  the  equation 
of  state.  A gas  constant  consistent  with  fully  dissociated  gaseous  Teflon  (i.  e. , 

2C  + 4F)  having  a mocular  weight  of  16.67g/mole  is  used  throughout. 

The  boundary  condition  on  magnetic  flux  density  at  the  ablation  surface  is  related 
to  the  condition  at  the  leading  edge  of  the  accelerating  plasma  at  each  time  step  for 
each  value  of  y through  Amperes  circuital  law 

where  integration  is  carried  out  over  any  closed  path  and  I is  the  current  flowing 
through  the  enclosed  surface.  Assuming  the  width  of  the  plasma  in  the  z-direction 
is  w,  and  denoting  the  position  of  the  leading  edge  of  the  plasma  at  each  value  of  y 
by  Xo(y)  the  integration  around  a rectangular  path  in  the  xz-plane  gives 

w -W&A  ♦ >*■)  -ty  V 


which  must  be  satisfied  at  each  value  of  y.  Bj  and  Bg  are  the  x and  z components  of 
magnetic  flux;  respectively.  No  variation  of  Bg  in  the  Z-direction  is  implied  in  the 
above  equation,  but,  nevertheless,  the  problem  does  assume  a certain  degree  of 
three  dimensionality  since  magnetostatic  field  components  will  and  do  exist  in  the 
space  around  the  plasma  because  the  acceleration  discharge  path  is  such  that  external 
fields  are  not  cancelled.  The  external  fields  are  due  to  the  flow  of  current  through 
the  plasma  itself  and  can  be  determined  once  the  current  density  throughout  the  plasma 
is  known.  The  relationship  which  must  be  evaluated  is 
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(3) 


The  values  of  B3  used  in  Equation  (2)  will  be  assumed  to  be  those  values  at  the 

center  of  the  arc  (z  = o).  Thus,  for  x = x . 

o 


where  f ^ +■  C*|- 

The  values  of  Bj  in  Equation  (2)  are  evaluated  at  the  outer  extrei  ities  of  the  arc. 
These  are  given  by 


Slice  B.osvS')- 

it  is  only  necessary  to  evaluate  one  of  the  above  and  integrate  it  from  x = o to  x = Xq, 
using  twice  the  value  of  the  resulting  integral  in  Equation  (2). 

Having  set  up  the  integrals  in  analytic  form  it  is  now  necessary  to  generate  a technique 
for  evaluating  them  numerically  at  each  time  step.  The  trapazoidal  rule  for  numerical 
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integration  is  used.  This  leads  to  rather  lengthy  formulas  and  for  brevity  only  the 
functional  form  of  a typical  integral  is  presented  below. 


5 y G^wxda -z 

The  above  procedure  leads  to  values  of  B3  at  x = Xo(y)  and  Bj^x,  y)  between  x - o 
and  x - xq  which  are  themselves  functions  of  B3(o,  y)  for  each  value  of  y.  This 
results  because  unless  B3(o,  y)  is  known  for  each  value  of  y it  is  not  possible  to 
calculate  the  current  density  at  the  point  x = Ax  using  a central  difference  approxima- 
tion. Thus,  if  M is  the  number  of  grid  lines  in  the  y-direction,  one  arrives  at  M 
simultaneous  linear  algebraic  equations  for  the  M values  of  B3  (o,  y);  y - - ^ MAy  to 

y = + £ MAy.  This  occurs  at  each  level  of  time  kAt.  The  numerical  solution  uses 
Gausian  elimination  to  solve  this  set  of  equations  for  the  values  of  B at  the  ablating 
wall.  Knowing  the  conditions  on  p , u,  v,  p and  B at  the  ablator  surface  allows 
calculation  of  the  dependent  variables  vector  W. 

The  leading  edge  of  the  accelerating  plasma  is  defined  as  that  point  at  which  the 
density  is  smaller  than  some  arbitrarily  small  value  which  is  held  constant  throughout 
calculations.  Ahead  of  this  point  all  values  of  p are  set  equal  to  zero  to  simulate  the 
true  vacuum  condition.  Similarly,  all  other  dependent  variables  are  zero  ahead  of  the 
plasma  except  of  course  the  magnetic  flux  density  which  must  decrease  in  magnitude 
from  its  value  at  the  leading  edge  to  zero  as  x approaches  infinity.  The  value  of  B 
at  the  leading  edge  is  obtained  at  the  same  time  the  values  at  the  ablation  surface 
are  obtained,  but  since  the  numerical  calculation  of  the  plasma  motion  in  the 
neighborhood  of  the  leading  edge  requires  three  additional  values  of  B,  i.  e. , namely 
those  at  Xq  + Ax  and  Xq  + 2 A x.  Equation  <4 ) is  utilized  twice  more  to  obtain  these 
required  values.  It  is  noted  that  in  truth  the  equation  governing  the  magnetic  flux 
density  ahead  of  the  plasma  is  wave-like  with  the  propogation  speed  equal  to  the 
speed  of  light.  Utilizing  a magnetostatic  field  equation  ahead  of  the  plasma  front  is, 
therefore,  only  an  approximation  to  the  real  state  of  affairs  but  is  an  excellent 
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approximation  when  the  front  velocity  is  not  relativistic,  as  is  the  case  in  this  device. 
Essentially  this  approximation  requires  that  the  unsteady  redistribution  of  flux  density 
ahead  of  the  plasma  occurs  at  essentially  infinite  speed  when  compared  to  the  speed 
at  which  the  plasma  leading  edge  is  moving  or  the  speed  at  which  the  value  of  B at  the 
leading  edge  is  changing.  This  is  in  fact  the  present  case  since  the  speed  of  light  is 
roughly  four  orders  of  magnitude  larger  than  the  plasma  velocity. 


Conditions  on  the  electrode  surfaces  are  considered  next.  The  fact  that  the  electrodes 
represent  impervious  surfaces  requires  that  the  normal  component  of  flow  velocity  be 
zero.  Regarding  the  tangential  component,  however,  a careful  evaluation  of  the 
continuum  boundary  condition  is  warranted.  Normally,  when  a fluid  is  assumed  to 
have  finite  viscosity  the  tangential  component  of  velocity  is  set  equal  to  zero  on  a 
boundary  surface.  It  is  known,  however,  that  under  certain  conditions  this  boundary 
condition  is  not  valid  and  a certain  amount  of  fluid  motion  tangential  to  the  surface 
does  exist.  When  these  conditions  exist  a so-called  "slip"  boundary  condition  is 
applied  as  opposed  to  the  usual  "no-slip"  continuum  situation.  The  parameter  which 
is  relevant  in  determining  whether  slip  or  no-slip  is  applied  is  the  local  mean  free 
path.  The  larger  the  mean  free  path,  the  less  valid  a no-slip  assumption  becomes. 
Velocity  at  the  surface  of  the  wall  is  given  by 


X 


where  X is  the  mean  free  path.  Under  conditions  of  low  wall  density,  the  velocity 
at  the  wall’s  surface  will  not  generally  be  small  as  compared  to  velocities  elsewhere. 


The  electrode  surface  temperature  will  be  considered  a constant  during  the  discharge 
period.  The  value  of  the  temperature  is  set  just  below  the  electrode  melting  point. 
For  copper  this  is  1335°K.  The  electrode  boiling  point  is  used  only  in  the  region 
which  the  plasma  exists.  Ahead  of  this  the  electrode  temperature  is  assumed  to  be 
room  temperature  (297°K).  The  boundary  condition  on  magnetic  flux  density  is 
obtained  from  conditions  which  must  be  satisfied  by  the  electric  field.  Namely,  if 
the  plasma  is  assumed  to  have  a relative  permitivity  of  unity,  the  tangential  com- 
ponent of  electric  field  at  the  electrode/plasma  interface  must  equal  zero.  This 
condition,  along  with  the  fact  that  v = o at  the  electrode  surface,  leads  to  the  result 
that  there  can  be  no  component  of  current  density  tangential  to  the  electrode  surface. 
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Hence,  the  normal  derivitivc  of  the  magnetic  flux  density  must  be  zero  at  the 
electrode  surfaces.  Utilizing  this  condition  together  with  internal  known  values  of 
H yields  the  value  of  U at  each  point  on  the  electrode  surfaces.  A scheme  similar 
to  that  used  at  the  ablation  surface  is  used  to  calculate  the  density  at  the  electrode 
surfaces. 

4.1.3  Plasma  Coupling  to  Discharge  Circuitry 

The  discharge  circuit,  consisting  of  the  energy  storage  capacitor  bank  and  strip  line 
assembly  connecting  the  bank  to  the  electrodes,  has  certain  capacitance,  resistance, 
and  inductance  associated  with  it.  The  plasma  represents  the  load  through  which  this 
circuit  is  completed  in  series.  Finite  plasma  conductivity  and  self  induced  magnetic 
field  yield  values  of  effective  plasma  resistance  and  inductance  which  must  be 
included  in  the  circuit  equations  so  as  to  simulate  the  physical  coupling  of  the  plasma 
load  and  capacitor  bank.  Both  the  plasma  resistance  and  inductance  are  functions  of 
time,  however,  since  the  plasma  motion  is  unsteady.  Moreover,  a nonlinear  coupling 
which  is  not  explicitly  obvious  in  either  the  circuit  equations  or  the  conservation 
equations  arises  because  the  boundary  conditions  on  mass  flow  rate  of  ablated  teflon 
and  magnetic  flux  density  are  functions  of  the  discharge  current  at  any  instant  of 
time,  which,  in  turn,  is  a function  of  the  temporal  variations  of  plasma  resistance 
and  inductance  up  until  that  time.  Thus,  it  is  necessary  to  solve  the  circuit  equations 
and  conservation  equations  simultaneously  to  include  this  physical  coupling  via 
boundary  conditions.  In  order  to  set  up  the  appropriate  circuit  equations  the  schematic 
diagram  shown  in  Figure  11  is  postulated  as  being  the  appropriate  equivalent  circuit. 

Expressions  for  the  plasma  resistance  and  inducatance  are  derived  from  the  basic 
definitions  of  these  quantities.  One  obtains  the  following  expressions  from  these 
definitions: 
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Figure  18 . Equivalent  Circuit  Schematic 

The  circuit  depicted  in  Figure  18  has  the  following  set  of  simultaneous  first  order 
differential  equations  associated  with  it: 

V(ft=  LS.+V+t|i+  <*L 

<**  (6  ) 

where  V (t)  is  the  potential  difference  across  the  capacitor  terminals.  Equations  ( 6 ) 
are  solved  using  the  second  order  accurate  Runge-Kutta  technique  with  the  integrals 
in  Equations  ( 5 ) evaluated  to  order  (Ax)2,  (Ay)2  using  the  trapezoidal  approximation. 

The  initial  conditions  on  Equations  (6  ) are  i = o and  (di/dt)  = Vq/Lq  at  t = o where 
L0  is  the  sum  of  Lc  and  I^(o)  determined  from  the  initial  arc  conditions  derived 
in  the  previous  section. 
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4.2  Governing  Equations 
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The  unsteady  fluid  dynamic  conservation  equations  with  two  dimensions  of  spatial 
freedom  may  Ik-  written  quite  generally  in  the  following  form: 

Wt  + fx  + S = FX  + Gy  + S (?) 

where  W is  a column  vector  containing  the  dependant  variables;  f and  g are  the 
advection  vectors  and  are  functions  of  W only;  F and  G are  the  dissipation  vectors 
which  are  functions  of  VV  and  its  first  spatial  derivitives;  and  S is  the  source  vector, 
a function  of  W only.  The  number  of  dependant  variables  (elements  in  W)  is  related 
to  the  particular  fluid  model  being  considered.  In  its  simplest  form,  that  is  for  a 
fully  ionized  one  species  plasma,  VV  will  contain  only  five  elements,  these  being  the 
density,  momentum  flux  in  both  directions,  total  energy  per  unit  volume,  and 
magnetic  flux.  Generalizing  the  model  to  a two-fluid  plasma  in  which  ions  and 
electrons  are  not  in  thermodynamic  equilibrium  but  are  moving  at  the  same  speed 
adds  an  additional  dependant  variable;  namely,  the  electron  energy  per  unit  volume. 

If  one  further  generalizes  the  model  to  include  a third  fluid  (neutral  gas  species) 
which  is  in  thermodynamic  equilibrium  with  the  ion  species  then  the  ionization 
kinetics  determine  whether  or  not  an  additional  dependant  variable  must  be  solved 
for.  Equilibrium  ionization  does  not  necessitate  addition  of  another  solution  variable 
since  the  chemical  composition  of  the  plasma  can  be  solved  for  using  the  algebraic 
equations  of  equilibrium  kinetics.  Essentially,  such  a kinetic  model  assumes  that 
all  ionization  reactions  are  occurring  infinitely  fast  in  both  the  forward  and  reverse 
directions  so  that  the  chemical  aggregation  of  the  plasma  can  adjust  to  changes  in  the 
thermodynamic  properties  at  each  point  in  space-time  as  quickly  as  the  flow  requires. 
The  assumption  of  finite  rate  ionization  does  require  that  an  additional  variable  be 
included.  Namely,  the  ion  or  electron  density  when  local  charge  neutrality  is 
assumed.  Moreover,  if  more  than  one  molecular  neutral  gas  species  is  present, 
and  finite  rate  dissociation  and  recombination  is  presumed  to  exist,  one  must  also 
include  the  neutral  gas  density  of  each  species  in  the  vector  W. 

The  model  used  herein  does  include  a neutral  gas  species  and  presumes  finite  rate 
ionization  but  does  not  include  the  existence  of  more  than  one  molecular  species. 

The  gaseous  products  produced  by  depolymerization  of  solid  tetrafluorethylene  (TFE) 
are  indeed  composed  of  many  different  probable  species  as  will  be  shown  in  the  next 
section  of  this  report,  but  these  various  neutral  molecules  have  been  more  or  less 
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lumped  together  for  reasons  which  are  explained  therein. 

The  resulting  set  of  eight  conservation  equations  (based  upon  the  above  discussion) 
are  defined  by  the  vectors  presented  in  Appendix  A of  this  report,  wherein  the 
effects  of  finite  electrical  and  thermal  conductivity  as  well  as  finite  viscosity  are 
manifested  in  F and  G whereas  finite  rate  ionization  and  therodynamic  nonequilibrium 
are  indicated  by  the  nonzero  elements  of  the  source  vector. 

4.3  Chemical  Model 

As  part  of  this  research  effort  a detailed  study  was  performed  in  which  the  probable 
chemical  kinetics  associated  with  gaseous  decomposed  TFE  was  analyzed.  The 
results  of  that  study  are  presented  in  this  section. 

With  regard  to  the  chemical  composition  of  the  plasma  produced  by  sublimation  of 
solid  TFE  very  little  information  is  available  in  the  literature.  Most  of  the  infor- 
mation which  is  available  has  been  generated  via  vacuum,  air,  or  inert  gas  pyrolisis 

11*10 

at  relatively  low  radiative  heat  flux.  Perhaps  the  most  relevant  piece  of  work 

17 

reported  in  the  literature  was  performed  by  Mathias  and  Miller  wherein  the  de- 
composition of  TFE  in  a glow  discharge  was  studied  using  helium,  helium  plus 
oxygen,  and  pure  oxygen  as  a carrier  gas.  It  was  determined  that  in  the  absence  of 
oxygen  the  major  constitutents  in  the  depoly mertzed  gaseous  products  are  CjF^ 
(85.22%),  C2F6  (6.43%),  n-C^Fg  (5.909?),  (2.08%),  and  CF4  (0.37%).  It  was 

also  determined  that  the  basic  sublimation  process  includes  unzipping  via  CF2 
radicals,  polymer  formation  via  C2 F4  excited  species,  and  CF2  radical  dispro- 
portionation to  account  for  the  F species.  The  reactions  thought  to  be  most  responsi- 
ble for  sustaining  the  discharge  are 

e + CF2  CF+  + F + 2e 

e + C2F2  C2F4  + 2e 

Because  of  the  low  ionization  potential  (10.3  ev),  C2F4  is  undoubtedly  a major  con- 
tributor of  additional  electrons  during  the  "avalanche"  phase  of  arc  development 
subsequent  to  breakdown. 

The  thermophysical  properties  relevant  to  ablation  of  TFE  have  been  the  subject  of 
more  thorough  investigation  than  has  the  chemical  kinetics  of  the  ablation  products 
(see,  e.g. , Refs.  16, 18).  These  studies  are  of  interest  since  they  have  provided 
some  relevant  information.  All  investigators  cited  are  In  agreement  that  the 


49 


FXUFtCMU-D 

f*L/BL.lC  « ()M»>4NV 


ablation  kinetics  are  first  order,  conforming  to  an  Arrhenius  formula  for  the  ablation 
rate  as  a function  of  temperature,  A general  picture  of  the  sublimation  process  con- 
sists of  solid  TFK  passing  through  a "gel"  state  at  327°C  in  which  the  polymer  bonds 
are  unzipped  (depolymerization)  and  the  monomer  is  converted  to  CF2  radicals.  This 

gel  remains  in  contact  with  the  solid  surface  and  becomes  gaseous  at  a temperature  of 

18 

640 °C  independent  of  the  rate  of  heating.  The  heat  of  ablation  quoted  in  the  litera- 
ture varies  from  333.4  and  424. 0 cal/g.  A mean  value  of  376.  8 cal/g  has  been  used 
herein.  Interestingly,  the  total  energy  required  to  ablate  the  . 148  mg  of  Teflon 
generated  during  a single  discharge  turns  out  to  be  only  about  one  quarter  of  a joule 
based  on  this  heat  of  ablation.  With  a discharge  energy  of  82J  this  is  an  insignificant 
fraction  of  the  energy  deposited  in  the  plasma. 

The  chemical  kinetic  model  postulated  herein  presupposes  that  the  products  of  ablation 
are  initially  composed  of  100%  monomer.  Justification  for  this  assumption  is  derived 
from  the  fact  that  even  during  the  extremely  low  energy  glow  discharge  of  Mathias  and 
Miller  only  about  fifteen  percent  of  the  products  were  observed  to  be  other  polymers. 
At  substantially  higher  heat  flux  Settlage  and  Siegle  observed  even  less  (3-5%)  amounts 
of  other  polymers.  It  is  therefore  more  than  likely  tht  all  polymers  are  decomposed 
at  the  level  of  heat  flux  experienced  by  the  ablation  surface  in  this  device.  The  gas 
composition  at  the  ablation  surface  is  therefore  presumed  to  be  100%  monomer  at 
640°C  and  is  invariant  with  time. 

Formulation  of  the  chemical  model  proceeds  by  writing  down  all  possible  dissociation 
reaction  leading  to  molecules  having  smaller  molecular  weight  than  CgF^.  Calcula- 
tion of  the  dissociation  energy  for  each  reaction  is  then  performed  by  taking  the 
difference  in  heats  of  formation  between  products  and  reactants.  Heats  of  formation 
for  the  molecules  considered  were  taken  from  the  JANAF  Thermochemical  Tables. 
These  are  referenced  to  zero  degrees  Kelvin.  The  result  of  the  above  procedure  is 
presented  below  where  possible  reactions  and  associated  energy  expressed  in 
kcal/mole  are  listed. 


(a) 
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C2  + 2F2  “ 350,5 

(b) 

11 

pH 

u* 

2CF2  - 94.0 

(c) 

o 

to 
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CF4  + C - 107.2 

<d) 

SF4^ 

CF3  + CF  - 107.9 

<e) 

C2F4^ 

2CF  + F2  - 301.4 
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(f) 
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Pm 
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u 

C2F2  + F2  ~ 104,2 

(g) 

CF 

C + F - 114.4 

(h) 

cf2  — 

CF  + F - 85.  3 

(i) 

CF3  ^ 

CF  + F2  - 193.5 

(3) 

cf3  — 

CF2  + F - 108.2 

(k) 

CF4  — 

CF3  + F - 115.1 

(1) 

o 

11 

CF2  + F2  " 186,5 

(m) 

CF4  — 

CF  + F2  + F - 308. 

(n) 

C2F2^ 

2CF  - 197.2 

<o) 

C2F2^ 

CF2  + C - 189.5 

Op) 

C2F2^ 

2C  + F2  - 389.2 

(q) 

F2  ^ 

2F  - 36.8 

On  the  basis  of  the  energy  requirements  for  the  above  reactions  it  appears  that  the 
existence  of  all  species,  i.e, , C,  F,  CF,  CF2,  CFg,  CF^  and  CgFg,  are  equally 
probable.  If  one  considers  the  overall  dissociation  energy  of  CgF^Si  2C  + 4F  a 
value  of  467.8  kcal/mole  is  obtained.  A molecular  weight  of  lOOg/mole  leads  to  a 
total  energy  requirement  of  4678  cal  to  completely  dissociate  one  gram  of  the 
monomer.  Hence,  a total  of  .692  cal  (2.9J)  is  required  to  completely  dissociate 
the  0.148  mg  of  ablated  TFE  generated  per  discharge.  It  would  therefore  appear 
that  a reasonable  approximation  to  the  chemical  kinetics  is  to  assume  that  the 
monomer  is  completely  dissociated  instantaneously  upon  being  ablated  because  of 
the  excessive  energy  being  deposited  into  the  resulting  gas.  This  approximation 
is,  however,  contingent  not  only  upon  the  energy  being  available  but  also  upon  the 
reaction  rate  being  sufficiently  high  with  respect  to  the  discharge  time  to  guarantee 
that  the  dissociation  will  proceed  to  completion.  Not  having  the  necessary  data  to 
compute  the  forward  rates  of  the  above  reactions  does  not  allow  verification  of  this 
approximation.  Moveover,  the  nonlinear  way  in  which  the  energy  is  deposited  into 
the  ablated  solid  leads  to  the  question  as  to  whether  or  not  sufficient  energy  per 
unit  mass  is  available  for  complete  dissociation  at  all  times  during  the  discharge. 
The  experimental  data  accumulated  throughout  this  study  suggests  that  after  the  first 
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cycle  of  the  discharge  oscillation  there  probably  is  not  sufficient  energy  being 
deposited  into  the  ablating  solid  to  completely  dissociate  the  gaseous  products  but 
the  high  density  and  temperature  near  the  surface  of  the  solid  propellant  and  the 
abundance  of  energy  available  for  dissociation  indicate  that  the  assumption  of  a fully 
dissociated  gas  is  probably  valid  at  least  for  the  first  cycle  of  the  discharge. 

Ionization  of  the  gas  mixture  2C  + 4F  is  considered  next.  The  first  four  ionization 
levels  for  fluorine  are  419. 6,  1226. 6,  2670.  9 and  4682. 7 kcal/mole.  For  carbon  gas 
these  are  430.2,  992.3,  2096.2,  and  3583.  1 kcal/mole.  Thermal  ionization  is 
postulated  as  being  the  major  mechanism  for  production  of  electrons.  The  corre- 
sponding rate  equation  for  this  process  is 


n.  = 

l 


= n.n  k,  - n k = -n 
i e f nr  n 


where  n.,  ng  and  nn  are  the  number  densities  of  ions,  electrons,  and  neutrals; 

respectively.  kf  is  the  forward  reaction  rate  constant  and  k is  the  reverse  reaction 
1 r 

rate  constant.  In  a static  equilibrium  situation  the  reverse  and  forward  rates  of 
reaction  must  be  equal  and  one  eliminates  having  to  know  one  of  the  rate  constants 
by  solving  for  it  in  terms  of  the  equilibrium  constant  kp.  The  following  is  then 
obtained: 


n. 
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= k. 


(n.n 
' x e 


- k n ) 
P n' 


where  kp  is  given  by  the  Saha  equation 


Kp  = 4.828  x 1015  G T/  exp  ( - ) 


where  G is  the  statistical  weight  (taken  to  be  1.5),  e,  is  the  ionization  potential,  T 
is  the  absolute  temperature  and  Kp  has  the  units  cm  . With  more  than  one  neutral 
and  ion  species  present,  the  above  rate  equation  is  written  for  each  different  species 
with  n.  and  nR  the  number  densities  of  the  ion  and  neutral  of  that  species  only 
whereas  ng  is  the  total  number  density  of  electrons.  Typically,  kf  will  obey  an 
Arrhenius  formula 


, oo  * 1 

kf  - CTe  exp  ( - -j^ — ) 

e 

where  C and  u>  are  constants.  For  neutrals  of  atomic  species  u:  usually  has  the 
value  ^ . The  value  of  the  constant  C is  only  approximately  determined  by  kinetic 
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theory  and  is  usually  adjusted  to  match  known  reaction  rates  by  multiplying  the  value 
obtained  from  kinetic  theory  by  a number  called  the  "steric  factor. " The  uncertainty 
in  the  value  of  C essentially  makes  it  an  unknown  unless  data  is  available  which 
allows  one  to  compute  the  proper  value.  No  such  data  is  available  for  ionization  of 
carbon  and  very  little  is  available  for  fluorine  ionization.  It  hardly  makes  sense, 
therefore,  to  consider  nonequilibrium  ionization  of  each  individually,  thereby  intro- 
ducing two  unknowns  when  one  can  approximate  the  true  state  of  affairs  by  assigning 
a single  forward  rate  to  electron  production  from  both  neutral  species  and  treat  the 
2C  + 4F  neutral  gas  mixture  as  a single  background  gas  having  an  average  molecular 
weight  of  16.67.  In  doing  so  the  problem  is  reduced  to  one  of  determining  a single 
value  of  C which  leads  to  the  correct  result.  The  implicit  assumption  is  that  the 
mole  fractions  of  C and  F remain  the  same  throughout  the  calculation. 

4.4  Numerical  Integration  Scheme 

The  computational  technique  to  be  used  for  solution  of  Equation  (7)  is  a second  order 
accurate,  two-step  explicit  numerical  integration  scheme.  The  scheme  has  been 
applied  successfully  in  the  past  for  solution  of  conventional  gas  dynamics  problems 
including  discontinuous  flows*^,  flows  containing  finite  rate  chemical  reactions  , 
viscous,  heat  conducting  flows^*,  and  radiative  flows  In  all  cases,  this  numerical 
technique  has  demonstrated  extremely  high  accuracy  when  compared  to  available 
analytic  solutions.  The  precise  form  of  the  scheme  is  given  in  Appendix  B. 

Since,  to  the  best  of  our  knowledge,  this  scheme  has  not  been  used  previously  to 

solve  MHD  problems,  a linearized  stability  analysis  was  performed  to  determine  the 

conditions  on  At  and  Ax  or  Ay  required  for  stability  of  the  calculations.  The 

analysis  consists  of  linearizing  the  finite  difference  approximation  to  Equation  (3)  so 

that  a single  equation  relating  W at  t + At  to  W at  t is  obtained  and  then  substituting  a 

Fourier  component  solution  of  the  form  a(t)  exp  (ittx  + By)  for  W (x,  y,  t).  One  then 

obtains  a system  of  five  simultaneous  linear  algebraic  equations  for  a(t  + At)  in  terms 

of  aft).  The  matrix  relating  these  two  vectors  is  called  the  amplification  matrix,  and  for 
' ' 2 
Equation  (3)  this  matrix  is  a function  of  the  arguments  (At/Ax),  (At/ Ay),  CTAt/Ax  , 

and  CTAt/Ay2;  where  CT  is  any  of  the  transport  coefficients  (i.e. , 1/a  u , v,  or  K). 

To  determine  whether  or  not  the  scheme  is  stable  in  this  linear  sense,  one  searches 

for  values  of  these  quantities  which  lead  to  eigenvalues  of  the  amplification  matrix 

having  a magnitude  less  than  one  over  all  space.  For  brevity,  the  details  of  this 

procedure  are  not  presented  here  and  only  the  results  are  discussed. 
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The  finite  difference  scheme  is  linearly  stable  in  the  case  of  zero  dissipation  for  all 
values  of  At/ Ax  satisfying  the  relationship 
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where  the  argument  on  the  right  hand  side  is  evaluated  at  each  point  in  the  flow  field 
at  each  level  of  time  and  the  minimum  value  is  used  to  compute  all  points  at  the  next 
level  of  time. 

This  result  is  a logical  extension  of  the  results  obtained  when  solving  non-dissipative, 
two-dimensional  unsteady  ordinary  fluid  dynamics  problems.  In  that  case  the 
characteristic  wave  speed  c is  the  speed  of  sound.  In  a MHD  flow,  the  speed  c is  the 
root  mean  square  of  the  speed  of  sound  and  the  Alfven  speed  (B  /up  )2»  In  the  case 
of  finite  dissipation  the  results  are  also  similar  to  those  obtained  for  ordinary  dissi- 
pative  fluid  dynamics  problems.  It  is  found  that  for  values  C-p  At/Ax'1  greater  than 
or  equal  to  approximately  10  stability  is  obtained  with  At/ Ax  given  by  Equation  (8); 
and  as  the  value  of  this  parameter  is  made  smaller,  the  value  of  At/ Ax  required  for 
stability  becomes  smaller  than  that  given  by  Equation  (4).  A discussion  of  the  reasons 
for  this  behavior  is  given  in  Ref.  21  where  it  is  also  concluded  that  to  obtain  the 
required  spatial  resolution  of  the  flow  field  and  retain  second  order  accuracy  of  the 
timewise  solution,  the  magnitude  of  the  time-step  must  be  the  same  or  smaller 
than  the  smallest  characteristic  diffusion  time. 
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5.0  RESULTS  OF  CALCULATIONS 

The  first  calculations  which  were  performed  were  done  using  the  one-dimensional 
unsteady  model  with  finite  electrical  conductivity  the  only  dissipative  mechanism. 

For  reasons  unknown  at  the  time  those  calculations  became  unstable  almost  imme- 
diately. Considerable  effort  was  exerted  in  varying  the  initial  conditions  to  deter- 
mine whether  or  not  an  inconsistent  set  of  initial  conditions  was  the  cause  of  the 
instabilities.  This  proved  not  to  be  the  case.  After  due  consideration  of  this 
numerical  problem  it  was  determined  that  the  actual  reason  for  development  of 
instability  was  an  unfortunate  choice  of  models.  Namely,  neglecting  the  displace- 
ment current  term  in  Maxwell's  equations  had  led  to  a set  of  MHD  equations  which 
were  physically  unrealistic  in  certain  regions  of  the  flow  field.  Precisely  how  the 
eroneous  physical  model  and  instability  are  related  is  discussed  below. 

One  of  the  approximations  to  Maxwell's  equations  which  is  normally  made  in  treating 
MHD  flows  is  to  neglect  displacement  currents  with  respect  to  conduction  currents. 
This  approximation  is  usually  a good  one  except  at  very  low  frequencies.  Disregard- 
ing displacement  currents  has  the  secondary  effect  of  eliminating  pure  electromagnetic 
wave  propagation  which,  in  most  flow  problems  of  interest,  is  not  a significant  effect. 
It  turns  out,  however,  that  in  a flow  such  as  the  one  being  analyzed  during  this  pro- 
gram where  the  advancing  plasma  is  interfaced  with  a true  vacuum  on  its  boundaries, 
the  propagation  of  electromagnetic  waves  ahead  of  the  plasma  is  important.  At  non- 
relativistic  plasma  velocities  the  effects  of  these  waves  on  the  overall  solution  to  the 
problem  can  still  be  neglected  but  in  a problem  wherein  the  characteristic  plasma 
wave  speed  can  approach  the  speed  of  light  the  assumption  of  negligible  displacement 
currents  must  be  reevaluated.  In  the  problem  at  hand  this  is  the  case  because  the 
basic  mechanism  associated  with  plasma  acceleration  is  an  expansion  which  converts 
both  thermal  and  magnetic  energy  to  kinetic  energy.  At  the  leading  edge  of  this 
expansion  (i.e. , the  vacuum/plasma  contact  surface)  the  plasma  density  must 
approach  zero  whereas  the  magnetic  flux  density  may  be  finite.  Thus,  the  character- 
istic wave  speed  (i.e.,  the  Alfven  speed)  at  the  contact  surface  will  approach  infinity 
if  the  displacement  current  is  neglected.  What  must  physically  happen  is  shown  in 
Figure  19.  The  fast  plasma  waves  will  asymptote  the  principal  electromagnetic 
waves  and  secondary  electromagnetic  waves  will  penetrate  the  plasma  asymptoting 
the  slow  plasma  waves.  This  will  occur  in  a continuous  fashion  as  the  collision 
dominated  (high  density)  plasma  expands  to  a free  molecule  flow.  In  the  original 
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Figure  19.  Typical  Wave  Propogation  Diagram 
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formulation  of  this  problem  neglecting  displacement  currents  allowed  the  wave  speed 
to  become  infinite  rather  than  being  limited  to  the  speed  of  light  and  the  instabilities 
encountered  were  simply  a result  of  not  being  able  to  choose  a value  of  At /Ax  smaller 
than  was  necessary  to  accommodate  the  physically  unrealistic  infinite  wave  speed. 

Inclusion  of  the  displacement  current  in  the  Biot-Savart  Law  introduces  considerable 
additional  complexity  to  the  governing  equations  due  to  the  appearance  of  additional 
terms  in  the  Lorentz  force  (JxB)  and  the  power  dissipation  ( E • _J).  Elimination  of 
E as  a dependent  variable  via  Ohm's  Law  and  Faraday's  Law  further  complicates 
matters.  One  finds,  however,  that  certain  reasonable  approximations  can  be  made 
which  eliminate  the  majority  of  these  additional  terms.  To  illustrate  this  the  one- 
dimensional unsteady  equations  will  be  considered. 

The  induction  equation,  obtained  by  differentiating  Ohm's  Law  once  with  respect  to  X, 
and  replacing  9E/9x  by  -SB/St  becomes 


Bt*  I"®'*'  W - 4BU> 

c 

which  differs  from  the  original  model  by  the  last  term.  The  Lorentz  force  takes  the 
following  form  with  displacement  currents  included: 


Jx  B = JB  = - ± E?  - C0luB2)t  - i fc  (U2B2)X 
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Only  the  first  term  on  the  right  hand  side  is  obtained  when  displacement  currents 
are  not  included.  Finally,  the  power  dissipation  term  becomes; 

E.J  = Ej  s — i—  (BB  ) - - (uB2)  - — B2  - E^  + — - (BEJ 

- - _ 2 ' x;  u ' 'x  2u  t 2 t rtu  ' t'x 

rru  x 


whereas  only  the  first  three  terms  appear  when  displacement  currents  are  negligible. 
The  first  fact  that  becomes  obvious  upon  looking  at  the  additional  terms  is  that  the 
electric  field  appears  explicitly  as  an  added  dependent  variable  which  cannot  be 
eliminated.  Hence,  Faraday's  Law  in  differential  form  would  be  included  in  the  list 
of  equations  to  be  solved  to  close  the  system.  The  second  and  more  important 
observation  is  that  the  differential  equations  to  be  solved  are  second  order  in  time 
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whereas  the  original  (no  displacement  current)  equations  were  first  order  in  time. 
Fortunately,  it  is  possible  to  simplify  the  resulting  equations  when  the  speed  of  the 
plasma  is  nonrelativistic.  In  order  to  make  this  approximation  we  note  that  within 
the  plasma  variations  in  time  are  of  the  same  order  as  the  product  of  velocity  and 
variations  in  space.  Mathematically  this  is  written  as 

0 (Wt)  = 0 (uWx) 

and,  therefore,  0 (W.J  = 0 (u2W  ) 

tt  ' xx' 

for  any  variable  W.  As  a result  of  this  it  is  seen  that  terms  of  the  form  wtt/c2  may 
be  neglected  as  being  small  compared  to  Wxx  because  (u/c)  is  much  less  than  one. 

It  is  also  noted  that  the  third  term  in  the  Lorentz  force  expression  may  be  neglected 
in  comparison  to  the  first  term  for  the  same  reason.  The  term  proportional  to 
(qB/H't  is  negligible  by  comparison  to  the  remaining  terms  for  typical  values  of  or. 
Near  the  contact  surface,  however,  <7  can  be  relatively  small  but  the  term  would 
vanish  identically  there  since  the  wave  equation  for  E and  B will  be  automatically 
satisfied  as  one  approaches  the  vacuum.  Indeed,  this  is  the  limit  of  the  MHD  equa- 
tions including  displacement  currents  when  rr  approaches  zero.  Additional  terms  in 
the  power  dissipation  term  are  systematically  eliminated  using  the  same  arguments 

and  the  additional  assumption  that  within  the  plasma  the  energy  stored  in  the  electric 
2 

field  (i.e. , cqE  /2)  is  small  compared  to  the  total  plasma  energy.  Thus,  the 
momentum  and  energy  equations  which  are  valid  approximations  to  the  generalized 
MHD  equations  under  the  assumptions  of  nonrelativistic  flow  velocity,  negligible 
electric  field  energy,  and  moderate  to  high  electrical  conductivity  are 

( ou  + f QuB2)t  = - (p*  + pu2)x 

E(.-  fr<E+P*>lx  <f  Bx) 

P x 

The  induction  and  continuity  equations  remain  unchanged  as  is  the  energy  equation. 
Note,  however,  that  the  momentum  equation  includes  an  additional  term  which  could 
not  be  neglected  on  the  basis  of  any  of  these  assumptions.  Rewirting  the  left  hand 
side  of  the  momentum  equation  as  follows: 
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it  is  immediately  obvious  that  the  additional  term  is  simply  proportional  to  the  square 
of  the  ratio  of  the  Alfven  speed  to  the  speed  of  light.  Under  conditions  wherein  the 
magnetic  flux  is  large  and  the  density  low  this  term  can  have  an  appreciable  effect  on 
the  solution  even  though  all  approximations  made  to  the  generalized  MHD  equations 
may  be  valid. 

The  procedure  used  in  calculations  performed  subsequent  to  the  above  reformulation 
is  as  follows.  The  solution  region  is  divided  into  two  domains;  between  the  ablation 
surface  and  the  plasma  leading  edge  wherein  the  MHD  equations  (including  the 
modified  momentum  equation)  are  utilized,  and  between  the  leading  edge  and  infinity, 
wherein  the  equations  governing  electromagnetic  wave  propagation  are  utilized.  The 
location  of  the  plasma  leading  edge  not  being  a known  function  of  time  a priori,  one 
must  also  compute  its  location.  Based  on  the  previous  arguments  a reasonable 
criterion  which  establishes  the  location  of  the  plasma  leading  edge  for  all  practical 
purposes  is  the  condition 

.2 


B 


LE 


LE 


For  values  of  BT  „ which  are  typically  on  the  order  of  one  T,  the  above  establishes  a 

-*11  3 

lower  limit  on  the  leading  edge  density  of  2.8  x 10  kg/m  . This  value  is  clearly 
beyond  the  acceptable  limits  of  continuum  theory  and  one  must  necessarily  be  cautious 
concerning  results  in  the  neighborhood  of  the  leading  edge.  The  two  solution  regimes 
are  "matched"  at  the  leading  edge  by  the  boundary  condition  on  B as  determined  in 
Section  4.1.2. 

Calculations  performed  using  the  procedure  described  above  allowed  stable  compu- 
tation to  proceed  much  farther  forward  in  time  than  was  possible  originally.  The 
separation  of  the  initial  arc  and  generation  of  a secondary  standing  arc  were  actually 
observed  in  the  numerical  solutions.  Results  of  a typical  run  are  presented  in 
Figure  20  where  B and  J are  plotted  as  a function  of  distance  at  four  different  times 
during  the  first  quarter  cycle.  As  postulated  through  analysis  of  the  experimental 
data,  what  appears  to  be  separation  of  the  initial  arc  is  actually  a splitting  of  the  arc 
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Figure  20.  Computed  Magnetic  Flux  and  Current  Density  Distributions 
During  the  First  Quarter  Cycle 
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into  two  separate  arcs  carrying  current  in  opposite  directions.  This  occurs  at  the 
onset  of  reversal  in  the  gradient  of  B at  the  ablation  surface  and  has  no  relationship 
to  the  time  at  which  the  discharge  current  passes  through  zero.  In  fact,  current  is 
already  flowing  in  the  opposite  direction  at  the  ablation  surface  when  discharge 
current  reverses,  eliminating  the  necessity  for  a "restrike"  of  the  arc.  The  reversal 
in  the  gradient  of  B at  the  ablation  surface  is  the  result  of  two  concurrent  phenomena. 
First,  as  the  discharge  current  passes  through  its  maximum  value  and  begins  to 
decay  so  does  the  magnitude  of  B at  the  ablation  surface.  Second,  the  magnetic  field 
energy  being  carried  away  from  the  ablator  by  the  plasma  is  continually  decreasing 
in  the  vicinity  of  the  wall  because  of  the  continual  conversion  of  this  energy  into 
kinetic  energy  via  the  MHD  expansion  process.  It  is  apprent,  however,  that  the 
characteristic  time  associated  with  the  decrease  in  the  value  of  B at  the  ablation 
surface  is  larger  than  the  characteristic  flow  time  in  the  vicinity  of  the  surface. 
Hence,  B decays  more  rapidly  at  the  surface,  the  rate  of  decrease  being  dominated 
by  the  overall  discharge  circuit,  than  it  does  as  a result  of  energy  conversion  near 
the  surface  (which  is  dominated  by  the  local  flow  characteristics). 

Numerical  instability  is  encountered  soon  after  arc  splitting.  The  reason  for  the 
instability  is  directly  associated  with  the  degeneration  into  two  separate  arcs,  which 
has  the  effect  of  producing  a low  density  region  between  the  standing  and  accelerating 
arc.  The  low  density  region  between  the  two  coexistent  arcs  is  caused  by  amplifica- 
tion of  the  magnetic  flux  density  between  the  arcs  tending  to  push  the  standing  arc 
against  the  ablator  surface  and  push  the  accelerating  arc  away  from  the  standing  arc. 
This  configuration  of  forces  exists  until  the  magnetic  flux  at  the  surface  of  the 
ablator  becomes  negative  (approximately  the  same  time  the  discharge  current 
reverses),  at  which  time  the  magnetic  pressure  gradient  within  the  standing  arc 
becomes  favorable.  During  the  period  of  time  that  this  force  configuration  exists 
no  mass  can  enter  the  region  between  the  two  arcs  and  extremely  low  densities  exist 
between  them.  As  soon  as  the  Alfven  speed  between  the  arcs  becomes  larger  than  the 
speed  of  light  the  calculations  again  become  unstable.  Note  that  this  problem  was 
circumvented  at  the  leading  edge  by  limiting  the  upstream  density  to  some  small  but 
finite  value  and  computing  the  magnetic  flux  ahead  of  it  using  standard  electro- 
magnetic wave  equations.  With  a similar  region  embedded  within  the  computational 
flow  field  such  a methodology  is  not  as  simple.  Indeed,  with  undefined  boundaries  on 
both  sides  of  the  region  the  problem  is  quite  unmanageable.  Finally,  after  many 
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attempts  at  resolving  this  problem  the  authors  were  forced  to  concede  and  concentrate 
on  resolution  of  the  acceleration  mechanics  for  the  first  quarter  cycle  only. 

Uncertainties  concerning  appropriate  values  for  plasma  conductivity,  scaling  length 
and  plasma  density  made  the  process  of  obtaining  good  quantitative  agreement  between 
the  computed  flow  field  and  experimental  data  a tedious  one.  A process  of  trial  and 
error  coupled  with  analytical  reasoning  was  undertaken  during  approximately  the  last 
year  of  this  program  with  that  objective.  Those  efforts  were  only  moderately  success- 
ful because  of  the  number  of  parameters  whose  unknown  values  could  encompass  a 
considerable  range.  Hence,  although  an  excellent  qualitative  picture  of  the  accelera- 
tion mechanism  was  obtained  no  absolute  quantitative  agreement  was  generated. 
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The  difficulties  encountered  in  the  numerical  formulation  of  boundary  conditions  for 
the  analytic  model  did  not  allow  calculations  to  proceed  much  beyond  the  first  quarter 
cycle  of  the  capacitor  discharge.  Nevertheless,  an  excellent  qualitative  picture  of 
the  acceleration  up  until  that  time  was  obtained  from  the  calculations  and  verified  by 
the  experimental  data  taken.  Good  quantitative  agreement  with  the  experimental  data 
was  not  obtained  because  of  the  uncertainty  in  the  values  of  several  important  para- 
meters. Despite  these  deficiencies  in  the  theoretical  modeling  certain  conclusiions 
concerning  plasma  acceleration  in  ablative  arcs  can  be  made  on  the  basis  of  what  it 
was  possible  to  accomplish. 

The  majority  of  the  discharge  energy  is  deposited  into  the  plasma  in  the  vicinity  of 
the  ablation  surface.  Energy  is  consumed  there  for  the  most  part  in  dissociating, 
ionizing  and  heating  the  ablated  mass  and  in  creation  of  magnetic  field.  Plasma 
acceleration  is  then  accomplished  via  an  MHD  expansion  in  which  magnetic  and,  to 
some  degree,  thermal  energy  are  converted  to  kinetic  energy.  The  process  may  be 
visualized  as  being  divided  into  two  regions;  an  energy  deposition  region  in  the 
vicinity  of  the  ablator,  and  an  energy  conversion  region  which  is  the  MHD  expansion 
fan.  There  is  a continuous  transition  between  these  two  regions  until  arc  splitting 
occurs,  a process  which  is  the  result  of  the  distribution  of  forces  as  the  discharge 
current  begins  to  decay  from  its  peak  value.  When  this  occurs  a portion  of  the 
initial  arc  remains  at  the  surface  of  the  ablator  and  the  remainder  continues  to 
accelerate.  Mass  is  retained  within  the  portion  left  standing  at  the  ablator  until  the 
magnetic  flux  changes  direction  at  the  surface,  at  which  time  the  standing  arc  also 
begins  accelerating.  The  region  between  the  two  arcs  has  very  low  density  because 
the  force  distribution  does  not  allow  for  the  flow  of  plasma  therein.  Upon  reversal 
of  the  discharge  current  the  above  process  appears  to  be  duplicated  with  the  second 
plasmoid  being  accelerated  in  much  the  same  fashion  as  the  first.  As  successive 
arcs  are  formed  and  accelerated  there  is  less  energy  available  to  each  and  the 
ultimate  velocity  which  each  plasmoid  reaches  is  therefore  lower.  The  fact  that 
less  energy  is  available  to  each  arc  as  the  discharge  progresses  is  offset  somewhat 
by  the  fact  that  each  successive  plasmoid  contains  less  mass,  which  maintains  the 
energy/mass  ratio  of  each  at  a high  value.  This  particular  observation  points  out  a 
very  basic  difference  between  plasma  acceleration  of  ablated  mass  products  as 
opposed  to  acceleration  of  plasmas  produced  in  gaseos  injection  type  devices.  That 
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is,  in  gaseous  propellant  devices  the  mass  flow  rate  of  gas  into  the  acceleration 
channel  is  constant  (being  limited  to  sonic  velocity  at  the  injection  orifice),  whereas 
ablative  plasma  generation  occurs  at  a rate  determined  by  the  discharge  current. 

The  rate  of  mass  addition  to  the  discharge  channel  at  any  given  time  is  therefore 
more  consistent  with  the  amount  of  energy  available  to  accelerate  that  mass  in 
ablative  plasma  accelerators  than  in  those  using  gaseous  injection  of  propellant. 

The  qualitative  understanding  of  the  acceleration  process  which  has  been  developed 
as  a result  of  this  program  does  not  immediately  yield  information  which  will  make 
it  possible  to  design  more  efficient  accelerators,  but  considerable  insight  into  the 
plasma  dynamics  of  such  devices  has  been  obtained  and  it  is  conceivable  that  further 
consideration  of  these  results  will  assist  in  future  development.  Plans  for  the 
immediate  future  are  directed  at  an  examination  of  the  plume  emanating  from 
millipound  thruster  systems  using  ablative  plasma  acceleration.  An  investigation 
of  the  plume  via  a pulsed  Langmuir  probe  to  be  developed  during  that  effort  will  be 
carried  out  and  a Monte  Carlo  stochastic  model  will  be  used  in  attempting  to  describe 
the  plume  analytically. 
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APPENDIX  A 

GOVERNING  EQUATIONS  - VECTOR  DEFINITIONS 


Nomenclature: 

f 

= 

Total  plasma  density 

P- 

= 

Density  of  species  «. 

u 

= 

x- component  of  velocity 

V 

= 

y-component  of  velocity 

E* 

= 

Energy  density  of  species  cx. 

E 

= 

Total  energy  density 

B 

= 

Magnetic  flux  density 

Pc 

= 

Partial  pressure  of  species  ot 

K 

= 

p*  + B2/2 

V1 

= 

2(2ux  - vy)/3 

V2 

- 

u + V 

y x 

V3 

= 

2(2vy  - ux)/3 

= 

Viscosity  coefficient 

<X 

= 

Electrical  conductivity 

K 

= 

Thermal  conductivity  of  species 

e« 

= 

Specific  internal  energy  of  species  oc 

= 

Collision  frequency  between  <X  and  ^ 

h> 

= 

Ionization  potential  (JAg) 

= 

Temperature  of  species 

mc 

= 

Mass  of  species  ot 

“i 

S 

Ion  production  rate 

67 


Dependent  Variables  Vector: 
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w = 


Advection  Vectors: 


P 

Pl 

P* 

pv 

Ee 

Ei 

B 


p * 

pv 

py 

jO*  + p\£ 

put/ 

/3UV 

3 =• 

fVpv1 

u(Ee4-f>«*) 

v(Ee+«) 

u(& 

uB 

vB 

Dissipation  Vectors: 


0 

0 

0 

0 

nv, 

’TV* 

w 

G = 

W, 

KS«+ 

»«,+  5^66, 

^(uv,  + 

^(uV^v^-v^ 

^4, 

L 
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Source  Vector: 


0 

mxnA 

0 

0 


A4tfer<5> 

0 


Algebraic  Relationships: 


^ ~ ft 

& = J2?p 
*«  tru.fi 

<?«=  I^*(E«/6)  - if-  CBViyw) 
ei=  fe  /fii)  ~ jf 


uVv1 
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APPENDIX  B 

FINITE  DIFFERENCE  SCHEME 


- vc- 

- W*(  <U-  9»V)  ♦ 

♦ txx  LFn:>r>-  f„ci + «ca 

x s A = At 

A-  “ Ax  Vii 

V/(x+n6x>i£  + tnA<4)&4-kA£') 

Ch*  3c*« 

4 MC--  O *yCvn-G«V-t) 

»C  C^V^l  “ 8^^ 

^*^*»*A  ~ ^ ^^Hry***!*  ^£ky«t ^n,»^t| 

^ ^»yn  = Ww.yN^  *;  CW^)^=  £^(W^m+-  ^ 

Jyj(  ^V**i  “^V  + H»tym2r7H»*iv*>) 
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processes  in  pulaed  ablative  plasma  accelerators  using  a continuum  MHD  approach 
as  the  basis  for  analytic  modeling.  Concurrent  with  analytic  modeling,  an 
experimental  diagnostic  effort  was  ateo  performed  so  that  calculated  results 
could  be  compared  to  the  measured  parameters  fof^evaluation  of  the  applicability 
aod-  oaeuracy  of  the  model.  Diagnostic  measurements  pegfst— d included  magnetic 
flux  density  using  small  coil  probes,  current  distribution  using  a Rog jwski  coil 
plasma  front  location  and  velocity  using  time  of  arrival  probes,  high  speed 
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photography  using  a high  speed  framing  camera,  and  electron  number  density  and 
temperature  using  Langmuir  probes.  Data  was  recorded  on  a digital  memory 
oscilloscope. having  appropriate  frequency  response.:*  An  unsteady  two-dimensional 
model  of  the  pTasma  motion  ^including  the  dissipative  effects  of  finite  plasma 
conductivity,  viscosity,  and  thermal  conductance^was  postulated.  The  model 
included  direct  coupling  of  the  external  capacitor  discharge  circuit  to  the  plasm  i 
via  the  time  variant  plasma  resistance  and  inductance. ^ A second  order  accurate 
two-step  explicit  numerical  integration  scheme  was  used  to Isolve  the  resulting 
seven  simultaneous  parabolic  partial  differential  equations  and  a second  order 
Runge-Kutta  technique  was  used  to  solve  the  two  simultaneous  second  order  ordinary 
differential  equations  for  discharge  current  and  potential.  ^Ablated 
(depolymerized)  solid  propellant  appears  to  be  dissociated  and  ionized  very 
quickly  with  respect  to  the  duration  of  the  discharge. Acceleration  of  the 
resultant  plasma  is  dominated  by  expansion  of  the  high  oensity  gaseous  products 
near  the  surface  of  the  ablator,  transforming  magnetic  and  thermal  energy  into 
kinetic  energy.  This  process  appears  to  be  repeated  during  each  quarter  cycle  of 
the  oscillatory  discharge  as  separate  plasma  bunches  are  generated.  The  process 
whereby  these  separate  plasmoids  are  produced  has  been  termed  "arc  splitting"  and 
is  the  result  of  a force  distribution  within  the  arc  which  is  dominated  by 
Lorentz  forces  and  the  ratio  of  characteristic  flow  time  to  characteristic 
discharge  time.  The  arc ' splitting  process  was  observed  numerically,  although 
quantitative  agreement  with  the  experimental  data  was  only  margional  because  of 
the  number  of  unknown  independent  parameters.  The  condition  of  large  Alfven 
speed  near  the  leading  edge  of  the  expansion  required  a unique  modification  to 
the  MHD  equations  to  insure  stability  of  the  calculation.  Future  work  will  have 
the  objective  of  describing  the  plume  emanating  from  pulsed  ablative  accelerators 
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